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Preface 



To reconcile classical theory of gravity with quantum mechanics is one of the most chal- 
lenging problems in theoretical physics. Einstein's General Theory of Relativity, which 
supersedes the Newton's law of universal gravitation, known since 17**^ century, is a geo- 
metric theory perfectly describing gravitational interactions observed in the macroscopic 
world. On the other hand Quantum Mechanics is indispensable for a description of mi- 
croscopic physics. At the energy scales corresponding to the Planck length ipi, due to the 
uncertainty principle, both quantum effects and gravitational interactions must be taken 
into account enforcing the necessity of quantizing also the gravitational fielc|^. 

The efforts to quantize directly the theory of gravitation using perturbative expansion 
meet serious theoretical problems, if such a theory is to be viewed as a fundamental theory 
and not only an effective theory. General relativity without matter field is perturbatively 
non-renormalizable at two loops [H [2]. Newton's constant G, which plays a role of the 
coupling constant of gravity, has a dimension [G] = — 2 in mass units and thus, as argued 
by Heisenberg, makes the theory power-counting non-renormalizable. Modification of the 
theory by adding higher derivative terms to the action can make it renormalizable but 
spoils the unitarity or the spacetime diffeomorphism invariance. 

A way to deal with this problem is to go beyond a conventional quantum field theory. 
String Theory attempts to unify all known interactions, the strong, weak, electromagnetic 
and also gravitation, treating general relativity as an effective low-energy limit of some 
fundamental theory [H]. 

Another set of theories postulate that quantum theory of gravitation has to be treated 
non-perturbatively. Loop Quantum Gravity is an approach to quantum gravity, which 
introduces new ways of treating gravity at the Planck scale by implementing the Dirac's 
procedure of canonical quantization to general relativity [H O [6] . Spin Foams [7] may be 
viewed as path integral formulation of Loop Quantum Gravity thus gaining some similar- 
ities with the model of Gausal Dynamical Triangulations (CDT), which is the subject of 
this thesis. 

From Loop Quantum Gravity originates a reduced model called Loop Quantum Gos- 
mology [8]. In the simplest case of isotropic geometries the model describes for example 
evolution of the scale factor using the effective Hamiltonian. Although in Causal Dy- 



The Planck length ipi = y 1.6 x 10 m where c is the speed of light in vacuum, G is the 

gravitational constant, and h is the reduced Planck constant. 
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Preface 



namical Triangulations we look at the same observable, there is a significant difference 
between the two approaches. In CDT we analyze the full quantum geometry and at the 
end we integrate out all degrees of freedom except the scale factor rather than introducing 
a reduction. 

The class of lattice approaches to quantum gravity may be divided into two groups. 
The first is based on fixed triangulations, containing the Regge- Calculus with edge-lengths 
as the dynamical variables O [ini E]. The second is using random triangulations, and 
includes Causal Dynamical Triangulations with fixed link lengths. The approach of Dy- 
namical Triangulations started with two-dimensional Dynamical Triangulations, a model 
of discrete random surfaces which appeared in various aspects. Let us mention only a small 
number of them. Its history begins with a combinatorial approach studied by Tutte[T^. 
By 't Hooft it was put in the context of two-dimensional Random Matrix Models in the 
large A^-limit and Quantum Chromodynamics [13] . This allowed to make a connection 
with a two-dimensional Euclidean quantum gravity \14:\ [T5] , also coupled to matter fields 

117] . Although analytical tools proved to be very powerful, there was still a problem 
with the definition of time which appeared to have a different scaling than the space. This 
inspired Ambj0rn and Loll to impose the causality condition which led to a formulation 
of the two-dimensional Causal Dynamical Triangulations [18]. The continuum limit of 
this model predicted the same two-loop amplitude as in the two-dimensional pure gravity 
theory obtained in the proper-time gauge in the continuum formulation [19]. However, 
even in two dimensions incorporation of matter fields was possible only by the use of nu- 
merical methods. The three-dimensional generalization was also widely studied [201 EI]) 
such models reveal a phase diagram with two phases similar to phases A and C of the 
four-dimensional model. 

In this dissertation we discuss the four-dimensional model of Causal Dynamical Trian- 
gulations without matter fields. It is a mundane lattice approach to quantum gravity using 
only standard quantum field theory and piecewise linear manifolds as a regularization, re- 
fraining from invoking exotic ingredients or excessive fine-tuning. There are premises 
based on the Renormalization Group (RG) approach that the asymptotic safety scenario, 
first justified by Weinberg [22], postulating the existence of a non-Gaussian ultra-violet 
fixed point is realized [23l [Ml [251 126] • The presented model attempts to define a non- 
perturbative quantum field theory which has as its infrared limit general relativity and on 
the other hand a nontrivial ultraviolet limit, which is in the spirit of the renormalization 
group approach, although the tools used by CDT are distinct from RG techniques. This 
indicates that lattice gravity theory may play the same role for quantum gravity as lattice 
field theory plays for quantum field theory. 

Subject of the thesis. The purpose of this thesis is to present recent results obtained 
in the framework of four-dimensional model of Causal Dynamical Triangulations. In par- 
ticular we give answers to questions like : how does a background geometry emerge dy- 
namically, what does it correspond to and how to describe quantum fiuctuations around 
the average geometry. 
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Structure of the thesis. The thesis is organized as fohows. 

In Chapter [1] we introduce methods of four-dimensional Causal Dynamical Triangula- 
tions and derive the foundations needed for further numerical and analytical computations. 
At the end of the Introduction, the author presents a list of results, where he participated 
during the preparation of this thesis. 

The logic of issues presented in the next four Chapters corresponds to a structure 
from top to bottom: we begin with a discussion of global macroscopic four-dimensional 
properties at long distances, and end with short-range fractal structure of slices. 

After reviewing in Chapter [2] the phase diagram of the theory we demonstrate the 
relation with Hofava-Lifshitz gravity and the physical importance of the de Sitter space. 

In Chapter [3] we prove that a background geometry which emerges dynamically cor- 
responds to the maximally symmetric solution of the minisuperspace model, namely the 
Euclidean de Sitter space. We study in detail the emerged geometry and show that it is 
genuinely four-dimensional and that in terms of lattice spacing it resembles an elongated 
spheroid. 

The work presented in Chapter S] aims to reconstruct the effective action describing 
quantum fluctuations of the scale factor around the semiclassical average. The resulted 
action is shown to agree with the discretization of the minisuperspace action. 

Chapter [5] contains results of measurements of Hausdorff dimension and spectral di- 
mension limited to hypersurfaces of constant time. Here we give a direct evidence for a 
fractal geometry of spatial slices. 

In the last Chapter [6] we describe in detail the Monte Carlo algorithm, used to obtain 
results presented in this work, and its implementation. 

Finally, in Conclusions we briefly discuss the main features of the four-dimensional 
model of Causal Dynamical Triangulations. 



Chapter 1 



Introduction to Causal Dynamical 
Triangulations 

The method of Causal Dynamical Triangulations (CDT) is a non-perturbative and back- 
ground independent approach to quantum theory of gravity. This model was proposed 
some years ago by J. Ambj0rn, J. Jurkiewicz and R. Loll with the aim of defining a lattice 
formulation of quantum gravity from first principles \28\ \30\ [29} [33] . The foundation of the 
model reviewed in this dissertation is the formalism of path-integrals applied to quantize 
a theory of gravitation. The standard way to define the quantum mechanical Feynman's 
path integral consist in introducing a discreteness of the time coordinate tj = e ■ i (for 
i = 1,...,N, e = ^) along non-classical trajectories of a particle which allows to ex- 
press the path integral by means of N ordinary integrals over particle positions and finally 
taking the continuum limit N ^ oo [341 [35] . 

The Causal Dynamical Triangulations method is a natural generalization of this dis- 
cretization procedure to higher dimensions. In the path integral formulation of quantum 
gravity, the role of a particle trajectory plays the geometry of four-dimensional spacetime. 
Causal Dynamical Triangulations provide an explicit recipe of calculating the path inte- 
gral and specify the class of virtual geometries which should be superimposed in the path 
integral. K. Wilson highlighted the significance of lattice field theory as an underlying 
non-perturbative definition of continuum quantum field theory [36]. Following this route, 
we hope that the lattice technique using causal dynamical triangulations as a regulariza- 
tion has the potential to play the same role in quantum gravity. Let us emphasize that 
no ad hoc discreetness of spacetime is assumed from the outset, and the discretization 
appears only as a regularization, which is intended to be removed in the continuum limit. 
The presented approach, has the virtue that it allows quantum gravity to be relatively 
easily represented and studied by computer simulations. 

The action. Classical theory of gravitation, General Relativity, in contrast with other 
known interactions describes the dynamics of spacetime geometry where the considered 
degree of freedom is the metric field g^y{x). The nonvanishing curvature of the under- 
lying spacetime geometry is interpreted as a gravitational field. The starting point for 
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construction of the quantum theory of gravitation is the classical Einstein- Hilbert action 
({—,+,+,+} signature and sign convention as in [371 [38]) 



where G and A are respectively the Newton's gravitational constant and the cosmological 
constant, M is the spacetime manifold equipped with a pseudo-Riemannian metric g^y 
with Minkowskian signature {—,+,+,+} and R denotes the associated Ricci scalar cur- 
vature [391 140j . Throughout this thesis we will use the natural Planck units c = h = 1, 
whereas the dependence on Newton's gravitational constant G will be kept explicit. We 
shall consider only closed manifolds, or more specifically, assume that the topology of Ai 
is X S^, and thus we pass over the Gibbons-Hawking- York boundary term. 

Partition function. Path-integrals are one of the most important tools used for the 
quantization of classical field theories. The path integral or partition function of quantum 
gravity is defined as a formal integral over all spacetime geometries, also called histories. 



In this expression we should integrate over geometries [g], i.e. equivalence classes of 
spacetime metrics g with respect to the diffeomorphism group Diffj^ on M. When 
integrating over metrics, one needs to divide out the volume of Diffj^, 



The partition function ()1.2p has a hidden dependence on coupling constants G and A. 
Knowledge of the generating function Z of a quantum field theory, after inclusion of the 
source terms, allows to calculate all correlation functions, i.e. vacuum expectation values 
of products of field operators, and provides the complete information about the theory. 

Causality. The underlying assumption of CDT is causality condition, which as we 
shall see, will have a significant impact on desirable properties of the theory. In |4H B2] 
Teitelboim advocated the requirement of causality by demanding that only those histo- 
ries for which the final three- geometry lies wholly in the future of the initial one should 
contribute to the path integral ()1.2p . The approach of Causal Dynamical Triangulations 
originates from this doctrine. In a gravitational path integral with the correct, Lorentzian 
signature of spacetime one should sum over causal geometries only. As an implication of 
causality, we will consider only globally hyperbolic pseudo-Riemannian manifolds, which 
allow introducing a global proper-time foliation. The leaves of the foliation are spatial 
three-dimensional Cauchy surfaces S. Because topology changes of the spatial slices are 
often associated with causality violation, we forbid the topology of the leaves to alter in 
time. Fig 1 1 . 1 1 illustrates a triangulation with imposed foliation which violates the causality 
condition. The spacetime topology may be written as a product M = I x T,, where / 
denotes an interval. For simplicity, we chose the spatial slices to have a fixed topology 




(1.1) 




(1.2) 




(1.3) 
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Figure 1.1: A visualization of a two-dimensional triangulation with a light-cone structure 
and a branching point marked. In Causal Dynamical Triangulations spatial slices are not 
allowed to split, which prevents singularities of the time arrow. 

E = 5^, that of a three-sphere. Moreover, we establish periodic boundary conditions in 
the time direction. Therefore, we assume spacetime topology to be = S*^ x S''^, where 
corresponds to time and 5^ to space. Such selection of Ai obviates the discussion of 
boundary conditions for the Universe. 

1.1 Causal triangulations 

The functional integration in (jl.2p is somewhat formal. To make sense of the gravitational 
path integral, the Causal Dynamical Triangulations model uses a standard method of 
regularization, and replaces the path integral over geometries by a sum over a discrete set 
T of all causal triangulations T. In other words Causal Dynamical Triangulations serve 
as a regularization of smooth spacetime histories present in the formal path integral (jl.2p 
with piecewise linear manifolds. Let us now define in greater detail, what we understand 
as a causal triangulation. 

The building blocks of four dimensional Causal Dynamical Triangulations are four- 
simplices. A simplex is a generalization of a triangle, which itself is a two-dimensional 
simplex, to higher dimensions. An n-simplex is an n-dimensional polytope with n + 1 
vertices. Concerning a n-simplex, sub-simplices of dimension n — 1 are called faces, while 
n — 2 dimensional sub-simplices are called hinges or bones. Each four-dimensional simplex 
is composed of five vertices (0-dimensional simplices), all of them connected to each other. 
It consists of 10 links (1-dimensional simplices) and of 10 triangles. The boundary of a 
simplex is built of five tetrahedral faces (3-dimensional simplices). Each four-simplex is 
taken to be a subset of a four-dimensional Minkowski spacetime together with its inherent 
light-cone structure thus the metric inside every simplex is flat. Fig. 11.21 presents a 
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visualization of four-simplices together with a light-cone sketch. As will be explained 
later, there are two types of simplices. 

A n-dimensional simplicial manifold, with a given topology, is obtained by properly 
gluing pairwise n-simplices along common faces. The neighborhood of each vertex (i.e. 
set of simplices sharing this vertex) should be homeomorphic to a n-dimensional ball. In 
this thesis we only consider four-dimensional simplicial manifolds topologically isomorphic 
with A4 = X S"^. A simplicial manifold takes over a metric from simplices of which it is 
built and which are equipped with a flat Minkowski metric. A simplicial manifold armed 
with a metric is called a piecewise linear space. In general such n-dimensional complex 
can not be embedded in M" which signifies a nonvanishing curvature. The curvature is 
singular and localized only at hinges, which in the four dimensional case correspond to 
triangles. 

As a consequence of the causality requirement we consider only globally hyperbolic 
manifolds with a proper-time foliation structure. In the Causal Dynamical Triangulations 
approach the spatial leaves of the foliation are called slices and are enumerated by a 
discrete time coordinate i. At each integer proper-time step i, a spatial slice itself forms 
a closed three-dimensional piecewise linear manifold T^^^(i) with a fixed topology of a 
three-sphere and induced metric which has a Euclidean signature. It is represented by 
a triangulation of /S^, made up of equilateral tetrahedra with a side length > 0. 

Two successive slices, e.g. at time steps i = t and i = t + 1, given respectively by 
triangulations T^^\t) and T^^\t + 1), are connected with four-simplices. The simplices 
are joined to create a four-dimensional piecewise linear geometry. Such object takes a form 
of a four-dimensional slab with a topology of [0, 1] x and has T^^\t) and T^^\t + 1) 
as the three-dimensional boundaries. A set of slabs glued one after another builds the 
whole simplicial complex. Such connection of two consecutive slices, by interpolating 
the "space" between them with properly glued four-simplices, does not spoil the causal 
structure. Null rays originating from tetrahedra in a slice t and propagating through 
successive flat simplices with Minkowski metric, always finally reach a slice t + 1 directly 
(without crossing other slices). This is true even for very twisted slabs. The triangulation 
of the later slice wholly lies in the future of the earlier one. 

Each vertex of a four-dimensional triangulation is assigned a discrete time coordinate 
i corresponding to the slice it belongs to. Because each simplex connects two consecutive 
spatial slices and contains vertices lying in both of them, there are four kinds of simplices: 
• type {4,1} with four vertices lying in a spatial slice T^^\t) and one in the next slice 
T^'^^t + 1), • type {1,4} obtained by interchanging the role of T^'^\t) and T^^\t + 1), 
i.e. with one vertex in the earlier slice and four in the next slice. • type {3, 2} with three 
vertices lying in a spatial slice T^^\t) and two in the next slice T^^\t + 1), • type {2, 3}, 
analogically, defined by interchanging T^^^t) and T^^\t + 1). The tetrahedra of T^^\t) 
are bases of four-simplices of type {1,4} and {4,1}. Sometimes we will not distinguish 
between types {1,4} or {4, 1}, and treat them as a common type {4, 1}, the same concerns 
{2,3} and {3,2}. Fig. 11.21 illustrates four-simplices of type {4,1} and {3,2} connecting 
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{4,1} {3,2} 



Figure 1.2: A visualization of fundamental building blocks of four-dimensional Causal 
Dynamical Triangulations - four-simplices (solid line). The simplices join two successive 
slices t and t + 1, and are divided into two types: {4,1} and {3,2}. The simplices are 
equipped with the flat Minkowski metric imposing the light-cone structure (blue drawing) . 
The dotted line illustrates the cross section of a simplex by a three-dimensional hyperplane 
placed between the slices, at time t + 1/2. 



slices t and t + 1. The simplices, as a subset of the flat Minkowski space, are equipped 
with the light-cone structure. 

Similarly, due to the causal structure, we distinguish two types of edges. The space-like 
links connect two vertices in the same slice, they have length ag > 0. The time-like links 
connect two vertices in adjacent slices and have length at- There occur no other types of 
links. In Causal Dynamical Triangulations, the lengths and at are constant but not 
necessarily equal to each other. Two different triangulations, indeed, correspond to two 
inequivalent geometries. By contrast, in the original Regge's model [A3\ lllj. the lengths 
of individual edges may vary, however, the gauge freedom in the edge lengths may result 
in an over counting of some triangulations. Let us denote the asymmetry factor between 
the two lengths by a 

al = a- al (1.4) 

In the Lorentzian case a < 0. The volumes and angles of simplices are functions of 
as and at and differ for the two types {4,1} and {3,2}. Because no coordinates are 
introduced, the CDT model is manifestly diffeomorphism-invariant. Such a formulation 
involves only geometric invariants like lengths and angles. The exact relations are derived 
in the Appendix A. 

Besides the slices with integer time coordinate, which are built only of equilateral 
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tetrahedra, we may also introduce slices with lialf-integer index i = t + i^- A horizontal 
section of a {4, 1} simplex with a three-dimensional hyperplane at half of its height is a 
tetrahedron. Similarly, by cutting {3, 2} with a hyperplane, we get a triangular prism as a 
boundary surface. The dotted lines shown on Fig. 11.21 illustrate the horizontal section for 
both kinds of simplices. A half- integer slice constructed of the two types of solids is also 
a closed three-dimensional piecewise linear manifolds and has a topology of three-sphere 

1.2 The Regge action and the Wick rotation 

The Einstein-Hilbert action (jl.ip has a natural realization on piecewise linear manifolds 
called the Regge action. Hereafter, we will denote the number of fc-dimensional sub- 
simplices by N^- Similarly, let mean the number of simplices of type {4, 1}, and A''32 
the number of simplices of type {3, 2}. Because we distinguish only two types of simplices, 
they have to sum up to the total number of simplices, 

Ni = Nii+N-i2. 

The total physical four-volume of a simplicial manifold T G T, is given by 

= j d^xV|detc/| = iV4iVol^^'^^ + iVgaVol^^'^^ 

where Vol^^'^-'' is a volume of a {4, l}-simplex and Vol'^-^'^J' is a volume of a {3, 2}-simplex. 
Both quantities are purely geometric and are proportional to and depend on a. Simi- 
larly, it can be shown that the global curvature 

j d'^xyJ\deig\R 

can be expressed using the angle deficits which are localized at triangles, and is a linear 
function of total volumes and the total number of vertices Nq. Based on the 

above arguments, the Regge action, calculated for a causal triangulation T, can be written 
in a very simple form, 

S[T\ = -Ko No[T] + N^T] + A [N^ilT] - 6No[T]), (1.5) 

where Kq, K4 and A are bare coupling constants, and naively they are functions of G, A 
and at,as- K4, plays a similar role as a cosmological constant, it controls the total volume. 
Kq may be viewed as inverse of the gravitational coupling constant G. A is related to 
the asymmetry factor a between lengths time-like and spatial- like links. It is zero when 
at = as and does not occur in the Euclidean Dynamical Triangulations. A will play an 
important role as it will allow to observe new phases. The derivation of the Regge action 
and explicit expressions for the coupling constants is comprised in the Appendix A. 

Discrete partition function. Causal Dynamical Triangulations provide a regular- 
ization method of histories appearing in the formal gravitational path integral ()1.2p . The 
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integral is now discretized by replacing it with a sum over the set of all causal triangu- 
lations T weighted with the Regge action (jl.Sp . providing a meaningful definition of the 
partition function, 

j;J-e^5[T]. (1.6) 

TeT 

Symmetry factor. The fraction I/C7- is a symmetry factor, given by the order Cj- 
of the automorphism group of a triangulation T. It might be viewed as the remnant 
of the division by the volume of the diffeomorphism group Dif fj^ present in a formal 
gauge-fixed continuous expression (|1.2p for the partition function Z. The factor I/C7- 
introduces the measure on the set of geometries, in the same way as does (jl.3p in the 
continuous case. In fact, we do not even have a mathematical characterization which 
geometries should contribute to the path integral. Therefore, the measure in ()1.6p is 
an attempt to define the quantum theory of gravitation. The next assumption is that 
piecewise linear geometries appearing in CDT are a dense subset in the set of geometries 
relevant for the path integral. There is no straightforward method to calculate Cj- for a 
general triangulation. However, we can easily deal with the factor Cj- by replacing the 
sum over unlabeled triangulations T € T in (|1.6p with a sum over labeled triangulations 
T G T. To each vertex of a triangulation T is assigned a unique label (e.g. an element 
of a set {vi} = {1, . . . , A'o[T']}). A triangulation is completely defined by the adjacency 
relations of simplices and all sub-simplices in terms of the vertex labels, preserving topo- 
logical restrictions. Two triangulations defined by the same set of vertex labels {vi}, are 
isomorphic if there exists a bijective map (j) '■ {^i} — {vi}, which maps neighbors into 
neighbors. The factor Cj- is defined as the number of such maps - among all A'^q! permu- 
tations, they generate the same set of neighbors, i.e. the same labeled triangulation. The 
number of labeled triangulations isomorphic with an unlabeled triangulation T (we may 
dress it with labels {vi} in a random way and apply the above definition of isomorphism 
class) is denoted as M[T], in other words it is the number of different labelings of T- From 
the definition of Cj- and MIT] follows. 



CTM[r] = No[r] 



The partition function ()1.6p can be written as a sum over labeled causal triangulations, 

Z=y -^e'^[^ = y -^e^^^-^. (1.7) 

The division by the number of vertex labels permutations ensures that we effectively sum 
over inequivalent parametrizations. Also from the numerical point of view it is much easier 
to consider labeled triangulations. The computer algorithm, described in Chapter [6l works 
on labeled triangulations and takes the symmetry factor Ct into account automatically 

isaisj. 

Wick rotation So far, it was impossible to evaluate the four-dimensional partition 
function ()1.2p or even its discrete counterpart ()1.7p using purely analytical methods. We 
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may however resort to numerical methods, namely to Monte Carlo techniques which allow 
to calculate expectation values of observables defined on piecewise linear manifolds. The 
advantage of the Causal Dynamical Triangulations approach, is that for a fixed size of 
triangulations understood as the number of simplices A'4, the number of combinations is 
finite, which in general allows to apply numerical calculations. Nonetheless, this number 
grows exponentially with the size, and the critical exponent coincides with the critical value 
j^cnt q£ bare coupling constant K4. Because of the oscillatory behavior of the integrand 
()1.2p or ()1.7p . we are still led into problems in defining the path integral, also the mentioned 
numerical techniques are not useful. We may evade this problem by applying a trick 
called Wick rotation, which, roughly, is based on the analytical continuation of the time 
coordinate to imaginary values, and results in the change of the spacetime signature from 
Lorentzian to Euclidean and a substitution of the complex amplitudes by real probabilities. 



In Causal Dynamical Triangulations, due to the global proper-time foliation, the Wick 
rotation is well defined. It can be simply implemented by analytical continuation of the 
lengths of all time- like edges, at — )• iat, 



This procedure is possible, because we have a distinction between time-like and space-like 
links. The Regge action rotated to the Euclidean sector, after redefinition applied in (|1.8p . 
gEuc _ _j^gLor ^ j^g^g gxactly the same simple form as its original Lorentzian version (|1.5p . 
An exact derivation of the Wick rotated Regge action is to be found in the Appendix A. 

Unboundedness of the Euclidean action. It is formally easy to perform the 
Wick rotation of the continuous Einstein-Hilbert action (jl.ip . However, the corresponding 
Euclidean action is unbounded from below, which is a consequence of the wrong sign of 
the kinetic term for the conformal mode. There are configurations with arbitrarily large 
negative value of the action, thus the Wick rotated contributions (jl.Sp to the path integral 
make it ill-defined. Also in the discretized model some triangulations may have very large 
negative values of the Regge action, but still finite due to the UV lattice regularization. 
The problem of infinities is revived when taking the continuum limit. Fortunately, in 
the non- perturb at ive approaches, like CDT, the partition function emerges as a subtle 
interplay of the entropic nature of triangulations, determined by the measure independent 
of hare coupling constants, and the hare action. The entropy factor may suppress the 
unbounded contributions coming from the conformal factor. This is exactly what happens 
in CDT. Together with a convergence of the coupling constants to their critical values, if 
such a point exists, the entropic and action terms should be balanced, and one hopes to 
obtain the proper continuum behavior. 

Let us note the importance of the causality assumption. The Causal Dynamical Trian- 
gulations approach arose as a modification of Euclidean Dynamical Triangulations, where 
no global foliation was imposed. This resulted in several problems. As proved by the 




(1.8) 



2 2 
a^ = a - a^ 



a > 0. 



1.3 The author's contribution to the field 
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computer simulations [65], at each point of the spacetime, the Universe could branch off 
creating a baby Universe. Inclusion of such degenerate geometries to the path integral, 
does not allow for a well-defined continuum theory. Because the Euclidean simplicial man- 
ifolds are bereft of the light-cone structure, one does not know how to rotate back to the 
Lorentzian signature and recover causality in a full quantum theory. Therefore in CDT 
we ab initio deal with the Lorentzian simplicial spacetimes in four-dimensions and insist 
that only causal well-behaved geometries appear in the regularized path integral. 

Expectation values. As a consequence of the regularization procedure and Wick 
rotation to the Euclidean signature, the partition function (|1.2p is finally written as a real 
sum over the set of all causal triangulations T (or labeled triangulations T), 

Z=y ^e-^[^] = y ^_e-^[^. (1.9) 

TeT TeT 

We should keep in mind, that the Euclidean Regge action S[T], as well as the partition 
function Z depend on bare coupling constants Kq, K4 and A. With the partition function 
p.9|) is associated a probability distribution on the space of triangulations P[T] which 
defines the quantum expectation value 

TeT 

where 0[T] = 0[T] denotes some observable. The above partition function defines a statis- 
tical mechanical problem which is free of oscillations and may be tackled in an approximate 
manner using Monte Carlo methods. Equation ()1.9p is a starting point for computer simu- 
lations, which further allow to measure expectation values defined by ()1.10|) and to obtain 
physically relevant information. 



1.3 The author's contribution to the field 

First, let us briefly summarize earlier results obtained within the CDT framework in four- 
dimensions. In the publication [30], the phase structure of the model was examined and 
three phases were found and characterized, namely phases A, B and C. They are also 
reported in the Chapter [2j Phase C became especially interesting from physical point of 
view. It was indirectly shown that a semiclassical background is generated dynamically. 
This background geometry corresponds to a four-dimensional de Sitter spacetime [3D1 132] • 
For that reason phase C is called also a de Sitter phase. Among others the spectral 
dimension, Hausdorff dimension and scaling properties were measured [31] . The repeated 
measurements are described in the Section 13.31 

Further, the author would like to depict his contribution to the development of four- 
dimensional Causal Dynamical Triangulations. The research described in this dissertation 
may be divided into two parts. The first part consisted of developing a computer software 
package allowing to perform numerical calculations within the CDT framework. The 
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author has created a set of programming tools ahowing to carry out simulations and 
data analysis. The main component of the package is the code generating statistically 
independent spacetime configurations according to a probability distribution (jl.lOp using 
Monte Carlo techniques. In general, simulations enable to measure the expectation values 
defined by (jl.lOp . The algorithm is based on a program written in FORTRAN by Prof. 
Jerzy Jurkiewicz in collaboration with Prof. Jan Ambj0rn and Prof. Renate Loll [29]. 
It was critically examined and fully rewritten in C by the author. As a part of the 
implementation of the Monte Carlo method the author has introduced and customized 
efficient algorithms and data structures, which are the domain of computer science. The 
used solutions, improvements and code optimization led to important acceleration of the 
code (ca. 7 times). Details of the implementation of the Monte Carlo algorithm and of 
the simulation process will be given in the Chapter EJ All numerical results presented in 
this work have been obtained by the author using computer programs written by himself. 
Here the author would like to pay tribute to the tools he used to complete the research. 
A vast majority of simulations presented in this work was carried out on a cluster at the 
Institute of Physics of Jagiellonian University. Parallel computations were performed on 
up to 64 cores. Some measurements required in total many CPU years of simulations. 
The author used the GNU Compiler Collection (gcc) and the Intel C++ Compiler (ice) 
to compile his source code. Wolfram Mathematica and Gnutplot allowed him to perform 
data analysis and create plots printed in this thesis. Some of the pictures were obtained 
with the help of Open Graphics Library (OpenGL) and CorelDRAW. 

Before performing measurements, one has to propose procedures determining what 
and how to measure. A construction of suitable observables which give physically relevant 
information is a highly non-trivial task. For instance, in quantum gravity one cannot speak 
about the absolute position of a point on the superposition of spacetimes. Moreover, they 
required specific numerical tools necessary to their measurement. 

The second part consisted of analysis and physical interpretation of results. The key 
element introduced by the author was a procedure eliminating the time-translational free- 
dom, which allows to superimpose configurations in an unambiguous way. The method 
is based on a unique designation of the triangulation center and is described in detail 
in Section 13.11 Breaking of the time-translational symmetry opened the door to further, 
more detailed analysis of the de Sitter phase. It allowed to introduce a new observable, 
namely the spatial volume at specific absolute time coordinate. Previously one could deal 
only with relative time positions. The crucial achievement was to directly show that the 
background four-dimensional de Sitter spacetime geometry emerges dynamically. In Sec- 
tion [3T] the author proves this statement. Examining the distribution of spatial volume 
for a specific slice, the author determined which region of the configurations is dominated 
by discrete effects and which can be treated as continuous. As shown in Section [3.2| the 
background geometry corresponds to the classical solution of a reduced minisuperspace 
model. The implementation of the Monte Carlo algorithm written by the author, allowed 
to repeat the measurements, e.g. of the spectral dimension, on larger configurations with 
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higher accuracy. In Section 13.31 the author presents the measured values of Hausdorff and 
spectral dimension, and based on the scaling properties shows that indeed the background 
spacetime is a well-defined four-dimensional Universe. Furthermore the author modified 
the method of measuring the spectral dimension consisting in performing the diffusion pro- 
cess to determine the propagator of scalar fields defined on a triangulation. This method 
is more efficient and accurate than the previously used heat-bath method. The results of 
these studies will be presented in a separate publication. In addition, improved numeri- 
cal tools allowed to investigate phase transitions with increased accuracy and to note the 
correspondence between the phase diagram of CDT and Lifshitz models, which form the 
basis for so called Hofava-Lifshitz gravity. Details of the correspondence will be given in 
the Chapter H 

In Section [3.41 the author examined in greater detail the global geometry of spacetime 
in de Sitter phase. It was shown, that background geometry resembles a four-dimensional 
spheroid, prolate in time direction, whose elongation in terms of lattice spacings depends 
on values of coupling constants. The author tested that it is possible to reach a spher- 
ical shape. Even for an individual configuration the appropriate volume functions did 
not deviate substantially from the average, suggesting finite volume fluctuations around 
background geometry. Also almost no fractal structure was observed with respect to the 
four-dimensional deflnition of distance. 

Another consequence of the elimination of the translational mode and the existence 
of background geometry is the possibility to perform a semiclassical expansion of the 
spatial volume around the mean value. This allowed for the determination of the effective 
action describing fluctuations of the three-volume. It should be noted here that we did 
not make any reduction of degrees of freedom, as is the case of cosmological models and 
minisuperspace model, but we integrated out all degrees of freedom except the scale factor. 
Moreover, the presented semiclassical expansion is truly non-perturbative and takes into 
account both the entropy factor and the hare action. However, it turned out that the 
resulting effective action shows a remarkable consistency with minisuperspace model. The 
author describes these results in Chapter HI 

In Chapter [5] the author conducted an in-depth analysis of three-dimensional slices 
- the surfaces of constant time with respect to the imposed foliation. He repeated the 
measurements of Hausdorff dimension and spectral dimension. In particular, he directly 
demonstrated the fractal nature of these slices. 
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This chapter is based on the following publications: (i) J. Ambj0rn, A. Gdrlich, S. Jordan, 
J. Jurkiewicz and R. Loll, "CDT meets Hofava-Lifshitz gravity", Phys. Lett. B 690, ^15 
(2010); (a) J. Ambj0rn, A. Gdrlich, J. Jurkiewicz and R. Loll, "CDT - an Entropic 
Theory of Quantum Gravity", [arXiv: 1007. 2560]. 

The Causal Dynamical Triangulations model is described by the Regge action ()1.5p 
which depends on a set of three bare coupling constants Kq, A and K^. All three coupling 
constants are nonlinear functions of parameters appearing in the continuous Einstein- 
Hilbert action, namely G and A, and the asymmetry factor a = which is a regularization 
parameter. Explicit relations are derived in the Appendix A. Parameters have following 
interpretations: Kq is proportional to the inverse bare gravitational coupling constant 
G. A is related to the asymmetry factor between lengths of time-like links at and space- 
like links Us. It is zero when at = and thus is not present in Euclidean Dynamical 
Triangulations (EDT). Here, A plays an important role as it allows to observe new phases. 
Finally, K4 acts as a cosmological constant A, it controls the total volume. For simulation 
technical reasons it is preferable to keep the total four-volume fluctuating around some 
finite prescribed value during Monte Carlo simulations. Thus needs to be tuned to its 
critical value, and effectively does not appear as a coupling constant. The two remaining 
bare coupling constants Kq and A can be freely adjusted and depending on their values 
we observe three qualitatively different behaviors of a typical configuration. The phase 
structure was first qualitatively described in a comprehensive publication [30] where three 
phases were labeled A, B and C. The first real phase diagram obtained due to large- 
scale computer simulations was described in [50]. The phase diagram, based on Monte 
Carlo measurements, is presented in Fig. 12. li The solid lines denote observed phase 
transition points for configurations of size 80000 simplices, while the dotted lines represent 
an interpolation. 

In the remainder of this Section we describe the properties of the phases, and discuss 
the phase transitions. 
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Figure 2.1: A sketch of the phase diagram of the four-dimensional Causal Dynamical 
Triangulations. The phases correspond to regions on the bare coupling constant Kq — A 
plane. We observe three phases: a crumpled phase A, a branched polymer phase B and 
the most interesting a genuinely four-dimensional de Sitter phase C. 
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Phase A. On the left: Snapshot of spatial volume N(i) for a typical configuration of 
phase A. The picture on the right presents a two-dimensional visualization of a 
corresponding triangulation. The time axis is set horizontally. 
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Phase A. For large values of Kq (cf. Fig. I2.ip the universe disintegrates into uncorrelated 
irregular sequences of maxima and minima with time extent of few steps. As an example of 
a configuration in this phase, the spatial volume distribution N(i), defined as the number of 
tetrahedra in a spatial slice labeled by a discrete time index i, is shown in the above figure 
(on the left). The picture on the right illustrates a two-dimensional causal triangulation, 
embedded in M^, with the time axis being horizontal. In this visualization, slices consist of 
cyclically connected edges and their circumferences correspond to spatial volumes. When 
looking along the time direction, we observe a number of small universes. The geometry 
appears to be oscillating in the time direction. They can merge and split with the passing of 
the Monte Carlo time. These universes are connected by necks, corresponding to minima, 
i.e. spatial slices not much larger than the smallest possible, which consists of five 
tetrahedra glued together. In the computer algorithm we do not allow these necks to 
vanish such that the configuration becomes disconnected. This phase is related to so-called 
branched polymers phase present in Euclidean Dynamical Triangulations (EDT) [27] . No 
spatially nor time extended universe, like the Universe we see in reality, is observed and 
the phase A is regarded as non-physical. 

160000 
120000 
§ 80000 



40000 
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Phase B. Snapshot of a spatial volume N{i) for an individual configuration in phase B 
and a two-dimensional visualization of a corresponding triangulation. 

Phase B. For small values of A nearly all simplices are localized on one spatial slice. 
Although we have a large three-volume collected at one spatial hypersurface of a topology 
of a three-sphere S^, the corresponding slice has almost no spatial extent. In very few steps 
it is possible to get from any tetrahedron to any other tetrahedron - along a path joining 
centers of neighboring tetrahedra. This is possible because many vertices, belonging to 
the largest slice, have extremely high coordination number, counted as the number of 
tetrahedra sharing a given vertex. The Hausdorff dimension is very high, if not infinite. 
In the case of infinite Hausdorff dimension the universe has neither time extent nor spatial 
extent, there is no geometry in a traditional sense. This phase corresponds to so-called 
crumpled phase also known from EDT. Phase B is also regarded as non-physical. 
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Phase C. Snapshot of a spatial volume N{i) for an individual configuration in phase C 
and a two-dimensional visualization of a corresponding triangulation. A typical 
configuration is bell-shaped with well-defined spatial and time extent. 

Phase C. For larger values of A we observe the third, physically most interesting, phase. 
In this range of bare coupling constants, a typical configuration is bell-shaped and be- 
haves like a well defined four-dimensional manifold (cf. above figure). The measurements 
of the Hausdorff dimensions confirm that at large scales the universe is genuinely four- 
dimensional [30]. Most results presented in this paper were obtained for a total volume 
of iV4 = 362000 simplices, and for Kq = 2.2, A = 0.6, K^"^ = 0.9. This point is firmly 
placed in the third phase C (cf. Fig. 12. A typical configuration has a finite time extent 
and spatial extent which scales as expected for a four-dimensional object. The averaged 
distribution of a spatial volume coincides with the distribution of Euclidean de Sitter space 
and thus this phase is also called a de Sitter phase. 

2.1 Phase transitions 

In this Section we will try to determine the order of the phase transitions. So far, there is 
a strong numerical evidence that the transition between phases A and C is of first order, 
while between phases B and C may be either a first-order or a strong second-order. Since 
we have at our disposal only results obtained by computer simulations, we can not verify 
the transition order with absolute certainty. 

Numerical determination of the position and order of the phase transition line demands 
very large resources of computation time. This is caused by several factors. To precisely 
determine the point of transition one has to perform dense sampling in the space spanned 
by the coupling constants. Often the observation of the phenomenon of phase coexistence, 
characteristic for the first-order transition, requires a very accurate tuning of bare cou- 
pling constants. In addition, we have to deal with the critical slowdown. Monte Carlo 
algorithms, which work well deep inside the phases, cease to be effective when approaching 
a transition point. The acceptance of Monte Carlo moves drastically decreases and the 
same moves tend to be repeated. This significantly extends the autocorrelation time. It 
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Figure 2.2: Transition between phases C (smaller Kq) and A (larger Kq) at Kq = 4.711 
and A = 0.6 for A'4 = 120k. We observe that configurations jump between two regimes, 
which is a strong evidence of first order transition. Top: Density plot of the spatial volume 
as a function of Monte Carlo simulation time (horizontal axis) and slicing time i (vertical 
axis). Darker colors mean larger spatial volumes N{i). Middle: Order parameter Nq/N^, 
conjugate to Kq, as a function of Monte Carlo time. Bottom: Distribution of the values 
taken by the order parameter Nq/N^. Double-peak structure is present. 

is therefore important to have efficient algorithms, which often must be customized for 
a specific transition. Also finite-size effects affect the position of the transition. This 
factor is particularly important for results presented in this Section. One can still expect 
shifts in the location of the line as A'4 — )• 00. A characteristic feature of first-order phase 
transitions is suddenness of the transition and the occurrence of the hysteresis with the 
change of coupling constants values in simulations. The latter however makes it difficult 
to find a precise location of a transition point. It should be noted that in numerical sim- 
ulations, which out of necessity work on finite configurations, one never observes a true 
phase transition which is related to singularities of derivatives of free energy . 

Let us first consider the A-C phase transition 0. As can be seen in Fig. 12.11 this transi- 
tion line can be approached by changing Kq, and keeping A, and for numerical reasons also 
A^4, fixed. A natural candidate for the order parameter is the variable conjugate to Kq, 

^The author would like to thank Samo Jordan for carrying out the measurements and providing the 
results. 
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namely the ratio Nq/N^. The graph in the middle of Fig. 12.21 shows the behavior of the 
order parameter N0/N4 as a function of the Monte Carlo simulation time at the transition 
line. As one can see, the order parameter jumps between two values which characterize 
different geometries in both phases. A density plot of the geometry shape, namely the 
spatial volume N(i), as a function of Monte Carlo time is shown at the top of Figure [2^21 
The quantity N{i) is plotted in the following way. On the horizontal axis is the simulation 
time r measured in multiples of Monte Carlo moves. On the vertical axis is the discrete 
proper-time i corresponding to the slice number. The value of A^^'^^ (i) is determined by the 
color of a point with coordinates (r, i). White color means a minimal three- volume, dark 
blue means the maximal volume and the intermediate values are properly graded. We 
can distinguish two types of behavior, either the three- volume gathers around a particular 
slice (blue marks) or it disintegrates over a whole time axis. There is a clear correlation 
between the type of geometry and the value of the order parameter. Large values of the 
ratio Nq/N^ correspond to a non- localized geometry, randomly distributed across the time 
axis. This is typical for configurations of phase A. Small values of the order parameter 
correspond to the situation where the geometry has a finite time extent and is localized 
around a single maximum. This is the expected behavior for the configuration of phase 
C. Clearly one can see that the nature of geometry changes abruptly, at fixed values of 
the coupling constants corresponding to the transition point. A centered configuration 
typical for the phase C, suddenly transforms into a disintegrated configuration typical for 
the phase A, after a while to reappear in another place. The graph at the bottom of 
Fig. 12.21 shows the distribution of values taken by the order parameter Nq/N4^. One sees 
two peaks, which correspond to different types of the geometry. This confirms the earlier 
hypothesis that configurations behave as if they were either in phase C or phase A and 
jump between them. Admittedly, the peaks are slightly blurred, but this is a consequence 
of finite-size corrections. With the increase of the total volume A4 — )• 00, peaks become 
sharper |50j . This suggests that, the A-C transition is of first-order. 

In a similar way we depict the measurements obtained for the B-C phase transition. 
Starting from phase C, we can approach phase B by varying a value of the coupling 
constant A (cf. Fig. 12. ip . Thus the variable conjugate to A (cf. Regge action (jl.Sp ). 
namely (A^4i — QNq)/N4, is used as the order parameter. Indeed, such choice of the order 
parameter, allows to observe similar effects as in the previous case. The graph in the 
middle of Fig. 12.31 shows the order parameter as a function of the Monte Carlo time. The 
same jumping behavior is observed as for the A-C transition. The graph at the top of 
Fig. 12.31 illustrates the geometry of configurations, by plotting the spatial volume N^'^\i) 
as a function of simulation time r. Again there is a strict correlation between the order 
parameter and the typical shape of the geometry. Large values of the parameter indicate 
that the system is in phase B. As can be seen from the top plot, in this case the universe 
has no time extent and the volume is wholly located at one slice. Small values of the 
parameter correspond to a system typical of phase C, which has a non-trivial time extent. 
The graph at the bottom of Fig. 12.31 plots the distribution of values taken by the order 
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Figure 2.3: Transition between phases C (larger A) and B (smaller A) at Kq = 2.2 and 
A = 0.022 for A'4 = 40A;. Although the configurations jump between two regimes, but 
the efi'ect gets weaker with increasing total volume A'4. Top: Density plot of the spatial 
volume as a function of Monte Carlo simulation time (horizontal axis) and slicing time i 
(vertical axis). Darker colors mean larger spatial volumes N{i). Middle: Order parameter 
(A'^41 — 6A''o)/A^4 , conjugate to A, as a function of Monte Carlo time. Bottom: Distribution 
of the values taken by the order parameter. Double-peak structure is present. 



parameter. As can be seen, the peaks of the distribution are even sharper and jumps are 
more evident. However, the strength of the first-order transition signal is suppressed with 
the increasing total volume A'4. Peaks become blurred and start to merge. Therefore, this 
result might be an artifact of too small configuration sizes. Moreover, the position of the 
transition depends on the total volume. If the system is near the B-C phase transition 
line, it happens that for small universe sizes the system is still in phase C, with time 
extended configurations. But, when increasing A'4, while keeping the coupling constants 
Kq and A fixed, the configuration width is decreasing, opposite to the expected scaling, 
and the configurations change their character and collapse to degenerated shapes of phase 
B. As already mentioned, due to the critical slowdown, we are not able to perform reliable 
measurements for larger volumes. So far, simulations did not answer the question whether 
the B-C phase transition is of first order, or if we observe a strong second-order transition. 

Because phases A and B are not regarded as physically relevant, we shall not study 
the phase transition between them. 
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2.2 Relation to Hofava-Lifshitz gravity 

Recently Petr Hofava proposed field-theoretical approach to quantum gravity inspired 
by the theory of a Lifshitz scalar. The modification of General Relativity introduces 
anisotropy between space and time, which makes the quantum theory power-counting 
renormalizable and potentially ultra-violet complete. The basic features of the theory, 
named Hofava-Lifshitz gravity, are only briefly mentioned in this Section, for details we 
refer the reader to the original works |52^ [531 EU • 

In the publication [52] Hofava points out some similarities between Hofava-Lifshitz 
gravity and Causal Dynamical Triangulations. Namely, the spectral dimension manifests 
the same short-range as well as long-range behavior in both theories. In the CDT frame- 
work, the spectral dimension is defined and calculated in the Section [HTHl It occurs, that at 
short distances it is equal to 2, while at large distances it is 4 as expected for the classical 
four-dimensional Universe. Identical scale dependence was obtain by Hofava in his theory. 
Another formal similarity of the two theories is the assumption of global time foliation, 
which refrains from equal treatment of time and space. In the de Sitter phase of CDT the 
anisotropy withers and the symmetry is regained. 

There are however further analogies between the two models. The Hofava-Lifshitz 
theory might be viewed as a generalization of the Lifshitz scalar theory to gravity. The 
theory of a Lifshitz scalar was first proposed to describe the tricritical phenomena and 
has a phase diagram which consists of three phases meeting at the Lifshitz triple point. 
There is a striking resemblance between the Lifshitz scalar phase diagram and the four- 
dimensional CDT phase structure. The Landau free-energy-density functional for a d- 
dimensional effective Lifshitz theory is given by 

F[<P{x)] = 020' + a40^ + • • • + C2{da(l)f + d2{d^cPf + 62(9^0)' + . . . , 

where (j){x) is the order parameter and /3 = l...m, a = m + 1 . . . d \55\ 156] . Distinction 
between directions a and /? allows for an anisotropic behavior and spatially modulated 
phases. A sketch of the phase diagram for the Lifshitz scalar is shown in Fig. 12.41 It 
is possible to make a one-to-one correspondence between the phases of CDT model and 
Lifshitz model if one assumes a following identification. We can qualitatively characterize 
geometries appearing in the three phases of CDT, and relate them to the Lifshitz mean- 
field order parameter (j){x). Let us introduce a qualitative notion describing the average 
geometry, and denote it as {geometry) . In phase C, we observe a genuine four-dimensional 
background geometry (cf. Chapter [3|) . The emerged Universe agrees with the classical 
notion of four-dimensional geometry, such that (geometry) > 0. While it might be con- 
sidered as a four-sphere S^, we can even argue that (geometry) = const. As stated at 
the beginning of this Chapter, configurations present in phase B have neither spatial nor 
time extent, thereby they have no geometry in the traditional sense, and we can write 
(geometry) = 0. Finally, in phase A, the configurations appear to be oscillating in the 
time direction, and \dt(geometry)\ > 0. So far, we do not know how to define a more 
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Figure 2.4: On the top: Sketch of Lifshitz scalar phase diagram with three phases heli- 
coidal, paramagnetic and ferromagnetic and the triple point. From left to right: the plot 
of the order parameter (p{x) respectively for helicoidal, paramagnetic and ferromagnetic 
phase. 

quantitative measure of the average geometry which could be identified with the order 
parameter. 

For m > 0, anisotropy is present and the system reveals three phases [§: 
Helicoidal. For (^2 < we have a modulated phase called helicoidal phase. The order 



■^Here we use the magnetic analogy to name the Lifshitz phases, but the phase diagram is apphcable to 
variety of other systems. 



30 



Phase diagram 



parameter is oscillating in /3 directions. For m = 1, we can identify this direction 
with the time present in CDT. In this phase we have > 0. Therefore, the 

helicoidal phase can be identified with phase A of four-dimensional CDT. The order 
parameter is illustrated at the bottom left picture of Fig. 12.41 

Paramagnetic. For d2 > 0, 02 > we observe a paramagnetic phase with a vanishing 
order parameter, (j){x) = 0. The order parameter is depicted at the bottom middle 
plot of Fig. 12.41 The paramagnetic phase can be identified with phase B present in 
CDT. 

Ferromagnetic. For d2 > 0, 02 < 0, we have a ferromagnetic phase. The order parame- 
ter is constant but nonzero, \4>{x)\ > 0, and is shown at the right of Fig. 12.41 Hence, 
the ferromagnetic phase can be identified with phase C of four-dimensional CDT. 

To complete the picture, let us consider the nature of phase transition lines. Most often 
transition between the ferromagnetic and paramagnetic phases in the Lifshitz theory is of 
second order. In CDT it corresponds to the B-C phase transition. So far, Monte Carlo 
simulations did not settle whether it is a first-order or a strong second-order transition. For 
Lifshitz scalar with m = 1 it might happen that the transition between the ferromagnetic 
and helicoidal phases is of first order, which would agree with the A-C phase transition 
in CDT. 

We can summarize, that the phase diagram of the Lifshitz scalar and the measured 
phase diagram of four-dimensional Causal Dynamical Triangulations model are strikingly 
similar. 



Chapter 3 

The macroscopic de Sitter 
Universe 

This chapter is based on the following publications: (i) A. Gdrlich, "Background Geom- 
etry in 4D Causal Dynamical Triangulations" , Acta Phys. Pol. B 39, 3343 (2008); (ii) 
J. Ambj0rn, A. Gdrlich, J. Jurkiewicz and R. Loll, "Planckian Birth of the Quantum de 
Sitter Universe", Phys. Rev. Lett. 100, 091304 (2008); (Hi) J. Ambj0rn, A. Gdrlich, 
J. Jurkiewicz and R. Loll, "Geometry of the quantum universe", Phys. Lett. B 690, 
(2010). 

In this Chapter we pass over local degrees of freedom of the quantum geometry, and 
reduce the considerations to volumes of spatial slices. In Causal Dynamical Triangulations, 
the causality condition is ensured by imposing on configurations a global proper-time 
foliation and keeping the topology of the leaves fixed. Due to the discrete structure, the 
successive spatial slices, i.e. hypersurfaces of constant time, are labeled by a discrete time 
parameter i. The index i ranges from 1 to T, and because of time-periodic boundary 
conditions the time slice z = T -|- 1 is cyclically identified with the time slice i = 1. To 
each vertex of a triangulation is assigned an integer time coordinate. All slices are built of 
equilateral spatial tetrahedra, each being a base of one simplex of the type {1,4} and one 
of the type {4, 1}. By construction they are glued in the way to form a simplicial manifold 
of a topology of a three-sphere S^. The topology of spatial slices is not allowed to change 
in time. 

3.1 Spatial volume 

The spatial volume N{i), also called the three- volume, is defined as the number of tetra- 
hedra building a spatial slice i = 1 . . . T. Because each spatial tetrahedron is a face shared 
by two simplices of type {1,4} and {4,1}, the three-volume N{i) sums up to the total 
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volume A'^tot equal 

Ntot = Y.N{i) = -N4,. (3.1) 

i=l 

Spatial volume N{i) is an example of the simplest observable providing information about 
the large scale shape of the universe appearing in CDT path integral. An individual space- 
time history contributing to the partition function is not an observable in the same way 
as a trajectory of a particle in the quantum-mechanical path integral is not. However, 
it is perfectly legitimate to talk about the expectation value {N{i)) as well as about the 
fluctuations around the mean. The average is defined by equation (jl.lOp . Due to the lack 
of analytical tools allowing to deal with the path integral (|1.9p of four-dimensional Causal 
Dynamical Triangulations, we have to resort to numerical computations in order to cal- 
culate expectation values. The Monte Carlo simulations generate a sequence of spacetime 
geometries, more precisely causal triangulations T, according to the probability distribu- 
tion (|l.iup . Configurations are then used to calculate the average. The measurement of 
the expectation values of observables is explained in detail in Chapter [UJ In this Chapter 
we show that the background geometry, described by the average {N{i)), emerges dynam- 
ically. In the next Chapter we study quantum fluctuations of N{i) around its average. 

All results reported in this Chapter correspond to one particular point of the phase 
diagram firmly placed in phase C, and given by the following values of bare coupling 
constants: Kq = 2.2, A = 0.6, volume A'^41 = 160000 and the time-period T = 80. In this 
phase, the plot of N{i) for an individual configuration is bell-shaped with a well-outlined 
blob. Below we describe how a typical triangulation of de Sitter phase is built and how to 
properly average N{i). Fig. 13.11 shows the spatial volume N(i) of a typical configuration 
as a function of a discrete time i. For the range of discrete volumes A'4 under study, the 
Universe does not extend over the entire axis, but rather is localized in a region much 
shorter than T = 80 time slices. We can distinguish three qualitatively different parts of 
the configuration (marked on Fig. 13. ip : 

• The stalk. The volume of a spatial slice is bounded from below by the volume of 
the minimal structure consisting of five tetrahedra connected to each other. Such 
a structure is the smallest non-degenerate triangulation of S^. This kinematical 
constraint ensures that the triangulation remains a simplicial manifold, in which, for 
example, two 4-simplices are not allowed to share more than one common face, and 
that it forms a connected space. The stalk is a region which is not much larger than 
the minimal structure. In Fig. 13.11 it is visible for i = —40 • • • — 22 and i = 22 . . . 39. 
The plot at the left of Fig. 13.21 presents the probability distribution Pi{N) which 
describes the probability that the three- volume of spatial slice i is equal to N{i) = N. 
The distribution is the same along all slices in the stalk. Thus, the mean volume 
is constant in this region. Apparently, there are very large lattice artifacts. The 
distribution Pi{N) shown in Fig. 13.21 splits into three groups: red for N{i) = 5 -|- 3k, 
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Figure 3.1: The spatial volume N{i) of a randomly chosen typical configuration in phase 
C {Kq = 2.2, A = 0.6). The bottom axis corresponds to shifted time index i with fixed 
position of the center of volume. The original time coordinate is presented on the top axis. 
We can distinguish three qualitatively different parts of the configuration: • the stalk with 
small fluctuations around the minimal slice structure, • the tail which volume is larger 
than the minimal value, but still very strong discrete effects are present, and • the blob 
which volume may be considered as a continuous variable. In the background, the blue 
curve plots the average spatial volume {N{i)) for the above values of coupling constants. 



green for N{i) = 5 + 3A; + 1 and blue for N{i) = 5 + 3A; + 2, where k is a natural 
number. Such behavior is an internal feature of the measure imposed on the set 
of triangulations and does not depend on the used Monte Carlo algorithm. The 
smallest possible volume is most probable and corresponds to the minimal structure 
with 5 tetrahedra, the slices larger by three tetrahedra are obtained by adding a 
vertex into one of the tetrahedra. To conclude, the slices in the stalk are usually 
at most several times larger than 5 tetrahedra and are dominated by discretization 
artifacts. 

• The blob. For slices with average volume much larger than the cut-off size {{N{i)) ^ 
5), the probability distribution Pi{N) becomes continuous and is well described by a 
normal distribution centered at the expectation value {N{i)). The right plot of Fig. 
13.21 presents Pi{N) for a central slice. The blob is a range of slices which fulfill above 
conditions, and it is evident for i = — 14 . . . 14 in Fig 13.11 Inside a blob lattice arti- 



34 



The macroscopic de Sitter Universe 



Probability distribution for the stalk 




Probability distribution for the blob 
Gaussian fit 




Figure 3.2: On the left: a probability distribution Pi{N) in the stalk. On the right: a 
probability distribution in the blob 



Probability distribution for the transition region 




Probability distribution for the transition region 




Figure 3.3: A probability distribution in the transition region inside the tail (for i = 17). 
For small the distribution splits into 3 families (left chart). For large N the split 
disappears but the distribution is highly asymmetric (right chart). 



facts become irrelevant and a semiclassical picture of Gaussian volume fluctuations 
around the mean is very reliable. 

• The tail. The intermediate region between the stalk and the blob is called the tail. 
Here, the volume is larger than the minimal value, but still very strong discrete effects 
are visible. It is present in Fig. 13.11 for |i| = 15. ..21. Fig. 13.31 shows the volume 
distribution Pi{N) for a particular slice in the tail. As shown in the right plot, the 
probability distribution is highly asymmetric and depends on the slice number. The 
left plot presents the zoomed part of the right picture, for A^ = 1 . . . 100. As in the 
case of the stalk, for A^ < 100 the three groups are present, and can be explained by 
the same mechanism. For larger volumes, the distribution decays exponentially and 
the split becomes unclear. 

The Einstein-Hilbert action (jl.ip is invariant under time translations t ^ t + 5. As a 
consequence, the same holds for the Regge action (jl.Sp under discrete time translations. 
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Because configurations are periodic in time, a straightforward average {N{i)) is meaning- 
less, as it would give a uniform distribution with a constant value, independent of i. From 
Fig. 13. H it is clear that in phase C the time translation symmetry is spontaneously broken. 
In order to perform a meaningful average of the spatial volume {N{i)), this should be 
taken into account and an appropriate time shift is needed. This is done by the proper 
centering of configurations, i.e. by fixing the position of the center of volume icv- Be- 
cause of the time translation invariance and periodicity, the following centering procedure 
is introduced: 

• Given a triangulation with spatial volumes N{i), we define the center of volume 
position relative to a slice i as 

T/2-1 

CV(i)= (5-0.5) •iV(i 5), (3.2) 

5=-T/2 

where we used the addition modulo T defined by z © j = 1 + mod(i + j — 1,T) and 
mod(i,T) is the remainder of the division of i by T. Because of the periodicity in 
time, the operator © is used instead of +. 

• We find the value of i G {1, 2, . . . , T} which minimizes |Cy(z)|, and denote this by 
icv- Due to the periodicity, |Cy(i)| has two minima, namely icv and icv © T/2. 
We choose the value of icv which is closer to the largest slice imax {imax roughly 
localizes the hloh center). 

• The Regge action is invariant also under time reflection z — ?• T + 1 — i. Therefore, 
we want the real center of mass to lie as close to icv + 0.5 as possible, i.e. in the 
interval {icv-,icv + !)• This explains the term —0.5 in ()3.2p . 

• Finally, we shift the definition of time, by an integer number of time steps, and fix 
the position of the center of volume icv to 0, 

N{i) ^ N{icv®i), i G -r/2 + l,...,r/2. 

Fig. 13.11 presents the spatial volume N{i) for a particular configuration as a function of 
the shifted time i (red line). The data represent the effect of the centering procedure, and 
the shifted time i is labeled at the bottom axis, while the original time index is shown at 
the top axis. The blue curve shows the average volume {N{i)). Evidently, the average 
position of the center of volume is localized at 0.5. 

Having fixed the position of the center of volume along the time-direction for all con- 
tributing triangulations, by applying the above shift procedure, we can now perform su- 
perpositions of configurations and define the expectation value {N{i)) as a function of the 
discrete time i. The average is measured using Monte Carlo techniques, 

(^(^))^;|E^^'^W' (3.3) 

k=l 



36 The macroscopic de Sitter Universe 



9000 

8000 

7000 

6000 
^ 5000 
^ 4000 

3000 

2000 

1000 


-40 -30 -20 -10 10 20 30 40 

i 

Figure 3.4: Background geometry {N{i)): Monte Carlo measurements for fixed A'41 = 
160000, = 2.2, A = 0.6. The best fit ()3.4p yield indistinguishable curves at given plot 
resolution. The bars height indicate the average size of quantum fluctuations. 

where the brackets (...) mean averaging over the whole ensemble of causal triangulations 
weighted with the Regge action ()1.5p and the expectation value is approximated be a 
sum over K statistically independent Monte Carlo configurations. The result for kih 
realization, is denoted by a superscript (.)*'*^\ where k = 1 . . . K. The details of evaluating 
the sum p.3p are described in Chapter [6l Fig. 13.41 shows the average spatial volume 
{N{i)) (red line) measured at a point in the phase C, Kq = 2.2 and A = 0.6. The height 
of the blue boxes visible in the plot indicates the amplitude of spatial volume fluctuations 
for each i given by cjj = {N{i)'^) — {N{i))'^. Results obtained by simulations show that 
the average geometry, in the blob and tail region, is extremely well approximated by a 
formula 

N{i) = (iV(i)) = H ■ cos^{i/W), (3.4) 

where W is proportional to the time extent of the Universe and H denotes its maximal 
spatial volume. The flt H • cos^ (i/W) is also plotted in Fig. 13.41 with a thin black line, but 
it is indistinguishable from the empirical curve. The background geometry given by the 
solution ()3.4p is consistent with a geometry of a four-sphere and corresponds to the 
Euclidean de Sitter space, the maximally symmetric solution of classical Einstein equa- 
tions with a positive cosmological constant [471 146j . This is one of the most important 
results described in this dissertation. We presented the direct evidence, that the back- 
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ground geometry of a four-sphere emerges dynamically. We did not put by hand any 
background. Moreover, neglecting the stalk, which by construction has a nonzero volume, 
we spontaneously end up with the topology, although we started with Ai = x S^. 



3.2 The minisuperspace model 

The shape of the three- volume N{i) = H ■ cos"^ (i/W) emerges as a classical solution of 
the minisuperspace model. The minisuperspace model appears for example in quantum 
cosmological theories developed by Hartle and Hawking in their semiclassical evaluation 
of the wave function of the Universe [57]. This model assumes spatial homogeneity and 
isotropy, which means that all degrees of freedom except the three-volume (scale factor) 
are frozen. In CDT model we have the opposite situation, no degrees of freedom are 
excluded, instead we integrate out all of them but the scale factor. Nevertheless in both 
cases results demonstrate high similarity. Neglecting internal structure of time slices, let 
us introduce a spatially homogeneous and isotropic metric on a Euclidean spacetime with 
X topology, 

ds^ = dr2 + a2(r)dO|, (3.5) 

where a(r) is the scale factor depending on the proper time r and dO| denotes the line 
element on S^. For such metric the physical volume of a spatial slice for a given time r 
equals 

v{t) = j dn^y^detgyi = 2Tr^a{Tf. 

and the scalar curvature is given hy R = ■^{1 — a'^ — aa), where the dot denotes a derivative 
w.r.t. T. The Euclidean version of the Einstein-Hilbert action (jl.ip is given by |37^ [38] 

SeWM = j dr dl73 \/dit^(i? - 2A). (3.6) 

The above action calculated for the metric ()3.5p up to boundary terms takes the form 

27r2 (■ 

S[a] = — — / dr i-Qah^ - Qa + 2Aa^) , (3.7) 
167rG J 

and is called the minisuperspace action. 

The Euclidean Einstein-Hilbert action (|3.6|) suffers from the unboundedness of the 
conformal mode, which is caused by the wrong sign of the kinetic term, as is reflected in 
the standard minisuperspace action (j3.7p . The same wrong sign is present in CDT, the 
Regge action (jl.5p is also unbounded from below, but the regularization makes the action 
finite for each individual triangulation. There is however a strong evidence [58[ 159]. that 
after integrating out all degrees of freedom, except the scale factor, which means taking 
into account the non-perturbative measure, one obtains a positive kinetic term in ()3.7p . 
As we shall see in the next Chapter, this is exactly what happens in four-dimensional 
CDT. The non-perturbative path integral over causal triangulations takes into account 
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both the entropy factor, counting the number of configurations, and the bare action. As a 
consequence, the effective action for N{i) is equal to the minisuperspace action (|3.7p with 
an opposite sign and removed unboundedness. 

For sufficiently small values of G the curvature term, generating negative sign of the 
kinetic term, will overcome the entropy factor, and histories with large time oscillations 
of the scale factor will dominate the path integral. In the Appendix A we derive relations 
between bare coupling constants, which show that Kq is approximately proportional to ^. 
For large values of Kq we expect to observe triangulations with large oscillations in time 
direction, and this is exactly what happens in phase A. 

Nonetheless, the overall sign of the action does not affect the classical solution of 
equations of motion. For the Lagrangian L = —aa? — a + ^Aa^, the Euler-Lagrange 
equation reads 

—— = + 2aa + Aa^ -1 = 0. 

da dr da 

The classical trajectory, solving the above nonlinear differential equation, is given by 

a{T) = R cos (t/R), i?= (A/3)-^/2. (3.8) 

Turning back to the spatial volume variable, the minisuperspace action p.7p can be rewrit- 
ten as ^ 

S[v] = y dr (^^ + /37;V3 _ 3A„^ . ^ = 9(2^2)2/3^ (3 9) 

The classical trajectory (|3.8p of physical volume equals 

v{t) = 27r2i?3 cos^ . (3.10) 

The physical volume v{t) describes the maximally symmetric space for a positive cosmo- 
logical constant, namely the Euclidean de Sitter Universe or a geometry of a four-sphere 
with radius R. This result is in agreement with the relation (|3.4p for N{i) found in 
numerical simulations. The de Sitter space emerges dynamically as a background geometry 
in the CDT model. In the next Chapter we give a direct evidence, that quantum fluctua- 
tions of the spatial volume N{i) are also governed by the minisuperspace action ()3.9p . up 
to the overall sign. 



3.3 The four dimensional spacetime 

In this Section, we prove that the Universe coming out in the CDT model is genuinely 
four dimensional. To support this statement we check the scaling properties and measure 
the spectral dimension of the appearing triangulations. 

Scaling dimension. Up to now, the measurements were performed for only one value 
of the total volume A'^tot- Keeping the coupling constants Kq and A fixed, which naively 
means that the geometry of simplices is not changed, we measure the expectation value 
iV(i) for different total volumes A'tot- 
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Figure 3.5: Average scaled spatial volume n{t) for a variety of total volumes A'^tot calcu- 
lated for the scaling dimension dn = 4. Measured in Monte Carlo simulations for Kq = 2.2 
and A = 0.6. We omit the error bars not to obscure the picture. The thin black line plots 
the fit n{t) = ^ cos^ {t/B), where B = 0.69. 

For the scaling dimension dn time intervals should scale as N^^" and the volume- 
independent time coordinate t scales as a function of a discrete time i as 

t^M-i, At^N;^'^. (3.11) 

In the continuum limit Atot — ^ oo, the time step should vanish At — ?• 0. Each configuration 
is composed of a stalk and a bulk. The kinematical constraint limits the minimal volume of 
slices lying in the thin stalk. Hence, one should use the effective total volume and subtract 
simplices which, out of necessity, lie in the stalk. Therefore, in the scaling relations (|3.1ip 
and (j3.12p appears the effective volume Atot = ^A'41 — Ts, where s 10 is the average 
spatial volume in the stalk. This subtraction has a very small effect for large configurations 
and vanishes in the continuum limit, but might be relatively important for the smaller 
ones. 

To compare spatial volume distributions N{i) for geometries with different volumes 
Atot) we introduce the scaled three- volume n{t). It has to be normalized and compatible 
with the time scaling ()3.1ip . This is satisfied by the following definition, 

n{t)^N;,l^'^Nii), n{t) = {nit)). (3.12) 



40 



The macroscopic de Sitter Universe 




Figure 3.6: Error of the overlap of scaled spatial volumes shown in Fig. 13.51 for various 
values of the scaling dimension dn- The point corresponding to the best fit is marked with 
a rectangle. 



For very large A'tot the time interval At is close to zero and in the continuum limit the 
sum over discrete time steps can be replaced by an integral 

J dt... o ^At.... (3.13) 

i 

Definition (|3.12p of the scaled spatial volume n{t) and equation ()3.ip for the total discrete 
volume A'tot provide the normalization condition 

/ n(t)dt = -^^iV(i) = l. 

Now it is possible to directly compare n{t) for various total volumes and check for 
which value of the scaling dimension du the overlap is the best. To estimate the optimum 
value of dn, we use a similar method as in [30] . For given values of coupling constants 
Kq, a and total volume Ntd, the average N{i) is measured and scaled in accordance with 
(IXTTTl and (fXT2D for fixed dn, 

(i, N{i)) ^ (^t = N;J^i, nit) = N~l^^Ni^)^ . 

The new series of points is then interpolated. Fig. 13.51 shows the scaled three-volumes 
n{t) using dn = 4.0 for several values of total volumes A'^tot- In order to obtain the best 
fit of dH, the error function is defined as follows. For a given value of the scaled time t, 
there is one value of the interpolated curve n{t) per each total volume A^tot- The error 
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function is defined as the variance of these values averaged over a number (of order T) 
equally distributed points t in the bulk region. The scaling dimension dn is chosen as 
the value corresponding to the minimum of the error function , i.e. to the best overlap of 
scaled spatial volumes n{t). Fig. 13.61 shows the dependence of the error function on cIh- 
The minimum corresponds to 

d// = 3.98 ±0.10. (3.14) 

The error of determination of cIh was estimated using the Jackknife method [601 EI] • The 
expected value dn = 4 is very close to the measured result, and is well within the margin 
of error. The result gives strong evidence, that the Universe which emerges in Causal 
Dynamical Triangulations is genuinely four- dimensional. 

Since n{t) is normalized and is obtained by the scaling of N(i) which is given by 
equation (|3.4p . it is expressed by the formula 

n(t) = A cos3 , (3.15) 

where B depends only on the coupling constants Kq and A, but not on A'^tot- For Kq = 2.2 
and A = 0.6, the measured values is i? ~ 0.69. The curve (j3.15p with adjusted B is drawn 
with a black line in Fig. 13.51 the fit is remarkably good. 

From equations ()3.12p and (|3.15p and the scaling dimension dn = 4 results a following 
expression for the three- volume N(i) 

As expected for a four-dimensional spacetime, the time extent Tuniv of the blob, measured 
in units of time steps, scales as Tuniv ~ t^B • N^^^. The expression p.l6p specifies ()3.4p 
and is only valid in the extended part of the Universe where the spatial three- volumes are 
larger than the minimal cut-off size. 

Physical volume. Let us relate the discrete spatial volume N{i) with the physical 
volume v(t) of hypersurfaces of a constant time. The classical solution of the physical 
volume v{t) is given by the formula (j3.10p . while the average discrete volume N{i) by the 
formula (|3.16p . Up to some factors they are expressed by the same function. Henceforth, 
we make a key assumption that the average configuration described by N{i) in fact has a 
geometry of a four-sphere S*^ given by v{t). The physical total four- volume of a four-sphere 
with a radius R equals 

= / v{r)dT = —R\ (3.17) 
J-^R 3 

On the other hand, the four- volume of a triangulations is given by (|A.2p . 

= Vol^4'^}a^iV4i + Vol^3.2}a4^^2 ^ (vol^^'^> + eVol^3,2}^ ^4^^^ ^ C^a'^N.^u (3.18) 

where ^ = jf^ and C4 = 2 ( Vol''-'^'^-'' -|- ^Vol'^-'^'^M , which is interpreted as the average 
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four-volume shared by one spatial tetrahedron. Here, a = is the cut-off length, i.e. 
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the lattice constant. The continuum time t defined by (j3.1ip and the discrete time i are 
proportional to the proper time r (cf. (|3.5p ). 



T = ^ft = ^fAti, At = Nj/\ (3.19) 

Therefore, a slab between slices i and i + 1 has a proper-time extent At and a four- volume 

v{t)At = v{T)^tAt = C^a^N{i) = Nll^C^a''n{t). (3.20) 

The above equation is consistent with formulas (|3.17|) and ()3.18p which determine the 
total four- volume of the emerging de Sitter space with a radius R equal 

R=[-^) a. (3.21) 

In the next Chapter, we will use the above relation to derive the renormalized value of the 
cut-off a. Comparing formulas ()3.10p and ()3.15p it follows that the proper-time extent of 
the de Sitter Universe is ttR, while in terms of the time t it is equal vri?, hence 

, T R , ^ 

V9rt = - = ^. (3.22) 

Assuming such scaling relations between physical and discrete volume (cf. (|3.20p ). and 
between proper and discrete times (cf. (j3.19p ). we ensure that the empirically derived 
formulas ()3.15p or ()3.16p describe a Euclidean de Sitter Universe for all A'^tot- 

Spectral dimension. Another quantity revealing information about the geometry 
is related to the diffusion phenomena, namely it is the spectral dimension ds- On a d- 
dimensional Riemannian manifold with a metric g^^{x), the probability density /j(x, xq; cr) 
describes the probability of finding a diffusing particle at position x after some fictitious 
diffusion time a with an initial position at cr = fixed at xq. The evolution of /j(x, xq; cr) 
is controlled by the diffusion equation 

9^/9(x, xo; cr) = Agp{x, xq; a), (3.23) 

with the initial condition 

p(x,xo;a = 0) = , V / '(x-xo), (3.24) 
Vdet5'(x) 

where is the Laplace operator corresponding to g^y{x). The return probability describes 
the probability of finding a particle at the initial point after diffusion time a. The average 
return probability P{o'), supplying a global information about the geometries, is given by 



P{<t) = {yj d'^xVdet<7(x)p(x,x;c7: 



where V4 = J d'^xy^det ^(x) is the total spacetime volume and the average is also per- 
formed over the ensemble of geometries. For infinite flat manifolds the spectral dimension 
ds can be extracted from the return probability due to its definition, 

4 = -2'^^!^. (3.25) 
d log cr 
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Figure 3.7: The spectral dimension ds of the Universe as a function of diffusion time 
a, measured for Kq = 2.2, A = 0.6 and N4 ~ 368k. The blue curve plots the average 
measured spectral dimension, while the highlighted area represents the error bars. The 
best fit ds{cr) = 4.02 — is drawn with red line. 

For Euclidean flat manifold W^, the spectral and Hausdorff dimensions are equal to the 
topological dimension, ds = dn = d. For the four-sphere 5"^, the spectral dimension 
£^5 = 4 for short diffusion times, while for very large times, because of the finite volume, 
the zero mode of the Laplacian will dominate and, with the above definition, dg will tend 
to zero. 

Definition ()3.25p is particularly convenient because it is easy to perform numerical 
simulations which measure the return probability. In the CDT framework, the spacetime 
geometry is regularized by piecewise flat manifolds built of four-simplices. Let us recall, 
that after the Wick rotation spacetimes appearing in the model are Riemannian manifolds 
equipped with the positive-definite metric tensor. The diffusion process can be carried out 
by implementing the discretized version of the diffusion equation ()3.23p 



where A denotes the time step and the sum is evaluated over all simplices j adjacent to i. 
Here variables iq , i and j denote the labels of simplices and the diffusion process is running 
on the dual lattice, i.e. the probability flows from a simplex to its neighbors. Since each 
simplex has exactly five neighbors, it is convenient to set A = 4 and the diffusion equation 



p{i, io; cr + 1) - p{i, io; cr) = A ^ {p{j, io;a) - p{i, iq; a 



)) 
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reads 

p(i,io;<T + 1) = -^/)(j,io;cr). (3.26) 

For a fixed triangulation, consisting of A^4 simplices, tlie procedure of evaluating p{i, io; a) 
is as follows [301 EI] : 

• We pick an initial four-simplex io. Because we are interested in the bulk features, 
and the volume distribution as a function of time is not uniform, the diffusion always 
starts from a simplex lying in the central slice icv- We impose the following initial 
condition ([3:2i]) 

p{i,io](T = 0) = 5ii^. 

• We iterate the diffusion equation ()3.26p and calculate the probability density /3(i, io; a) 
for consecutive diffusion steps a. 

• Finally, we repeat the above operations for a number of random starting points iq 
adjacent to the central slice (because a simulation of the diffusion process is quite 
time-consuming it is repeated K = 100 times per configuration) and calculate the 
average return probability 



1 ^ 

Pier) = P(^o,«o;f^) 



io=l 

In numerical simulations the return probability P{(t) is averaged over a number of trian- 
gulations (~ 1000) and the spectral dimension ds is calculated from the definition ()3.25p . 
Fig. 13.71 shows the spectral dimension ds as a function of the diffusion time steps a , in 
the range 40 < a < 500. For small values of cr (< 30) lattice artifacts are very strong and 
the spectral dimension becomes irregular. Because of the finite volumes of configurations, 
for very large cj(^ 500), the spectral dimension ds falls down to zero. In the presented 
range, the measured spectral dimension ds is very well expressed by the formula 

dsia) = a- = 4.02 - (3.27) 

c + a 00 + a 

where variables a, b and c were obtained from the best fit. As observed, the spectral 
dimension depends on a diffusion time, and thus it is scale dependent. Small a, means that 
the diffusion process probes only the nearest vicinity of the initial point. Extrapolation of 
results gives the short- distance limit of the spectral dimension 

ds{a ^ 0) = 1.95 ±0.10. 

In the long-distance limit the spectral dimension tends to 



ds{a ^oo) =4.02 ±0.05. 
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The measurements presented here were performed on configurations twice as large as 
those used in j3H [30] . but results are the same. The short-range value of the spectral 
dimension ds = 2, much smaller than the scaling dimension du, suggests a fractal nature of 
geometries appearing in the path integral at short distances. This statement is supported 
by results presented in Chapter [5] where we give a direct evidence of the fractal nature of 
spatial slices. At long distances ds = 4, and configurations resemble a smooth manifold, as 
we show in the next Section. Amazingly such nontrivial scale dependence of the spectral 
dimension of the quantum spacetime, the same infra-red {ds = 4) and ultra-violet {ds = 2) 
behavior, is also present in the Hofava-Lifshitz gravity [52] and Renormalization Group 
approach [52] . 
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In the previous Sections we identified the background geometry with the geometry of a 
four-sphere. The purpose of this Section is to describe the geometry in more detail, with 
a particular emphasis on approaching the B-C transition line by varying the coupling 
constant A. We show that typical configurations in phase C resemble a prolate spheroid 
stretched in the time direction, and do not exhibit a fractal nature. The three-volume 
N{i) was defined as a number of tetrahedra in the slice i, of a centered triangulation. For 
convenience, here we will use a slightly different definition of volume. The spacetime slab 
between slices i and i + \ has a four- volume denoted by Ni{i), which includes simplices 
of all types {4,1} and {3,2}. The two definitions are approximately equivalent, since 
N4{i) oc N{i). Fig. 13.81 presents the average slab four-volume A^4(i) for three different 
values of the coupling constant A. In addition to A^4(«), we introduce the radial four- 
volume A^4(r). The shell of radius r is defined as a group of simplices whose distance from 
some initial simplex equals r. This distance is defined as the length of the shortest 
geodesic connecting the two points. Although we may determine the exact geodesic on a 
piecewise linear manifold, we use a discrete geodesic distance defined as the length of the 
shortest path between successive centers of neighboring four-simplices. We expect that 
our definition of a geodesic will in the continuum limit, i.e. on scales sufficiently large 
compared to the cut-off scale, lead to the same geometric results as for exact geodesies. 
The number of simplices contained in a shell of radius r is denoted as N4{r). The initial 
four-simplex is located at r = 0, thus A'4(r = 0) = 1. Each simplex has five neighbors: 
A^4(r = 1) = 5. Moving out we find successive shells. Moreover, we introduce a partition 
of simplices belonging to a slice r with respect to the time slab i they lie in. The number 
of simplices lying simultaneously in shell r and time slab i is denoted as N4{i, r) and the 
following relations are satisfied. 



iV4(i) = Yl r), Ni{r) = ^ N^, r), 
T i 
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Figure 3.8: The average time slab four-volume (A^4(i)) shown for A = 0.1,0.25,0.6 and 
Kq = 2.2. N/^^i) denotes the number of all four-simplices contained between slices i and 



summing up to the total volume 

iV4 = ^Ar4(i) = ^iV4(r). 

i r 

The fractal structure. To examine a fractal nature of a specific configuration, we 
define a diffusion tree, which measures the connectedness of radial shells. Two simplices 
in a shell of radius r are considered to be connected if one can find a path connecting them 
passing only through simplices from the outer shells r' > r. A diffusion tree is created in 
the following way: 

1. For a given configuration we choose the initial simplex iq from the central slab. 

2. We put this simplex to the shell r = and remove it from the triangulation. 

3. We find the first shell, i.e. all survived simplices which were connected to the original 
simplex. We remove all simplices in this shell from the configuration. 

4. We find the next shell, i.e. all survived simplices which were connected to at least 
one simplex from the previous shell. We remove this shell. 

5. We repeat step 4. until no simplices remain. After each removal, we check a number 
of disconnected parts of the remaining triangulation. 
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Figure 3.9: Graphs of the diffusion trees. From left to right the figures correspond to 
spacetime configurations measured for A = 0.10,0.25 and 0.60. Each vertex represents a 
connected group of simplices in one shell. The top node contains only the starting four- 
simplex. The distance from the top node equals the shell radius r. For A = 0.25 and 0.60 
the bifurcation of the branch at radius r^if is observed. 

6. A set of simplices belonging to the same shell and the same connected component 
(they may be connected through simplices in further shells, but not by the removed 
ones) is represented as a vertex on the graph. 

7. A vertex is linked to another vertex in the previous shell, if the corresponding com- 
ponents share neighboring simplices. Except for the original simplex, each vertex 
has exactly one parent. 

By construction such a graph forms a tree, and the procedure of creating it imitates 
the diffusion process, as the shells follow the front of diffusion. Fig. 13.91 illustrates the 
diffusion trees obtained for typical configurations corresponding to three points with A = 
0.10,0.25,0.60 and Kq = 2.2. Only vertices which volume, together with sub-leaves, is 
larger than 40 are visible. Each node represents a connected component of a shell and its 
distance from the top node is equal to the shell radius. Vertices are drawn as circles whose 
diameter depends on the size of the components they represent. 

The average volume {N^^i)) for different values of A is presented in Fig. 13.81 Similarly 
as N{i), shown in Fig. 13.41 Ni{i) is expressed by formula N4^{i) oc cos^ i/W. This behavior 
corresponds to a four-dimensional spheroid, with three spatial radii of equal length. Al- 
though we know the time extent of the spheroid (cf. Fig. 13. Sp . so far we can not say much 
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Figure 3.10: The contour plots of {N4,{i,r)) for A = 0.10,0.25,0.6 and Kq = 2.2. Darker 
colors mean larger values of {N4{i,r)). The horizontal axis denotes the discrete time i. 
The vertical axis denotes the radius r. 



about the length of the equator (or equivalently spatial radii) as we have no information 
about the geometric nature of spatial slices. An individual configuration usually is given 
by a deviation of the mean geometry. 

Studying the properties of diffusion trees, illustrated in Fig. 13.91 we indicate that 
typical configurations, corresponding to the two largest values of A = 0.25, 0.6, resemble 
a prolate spheroid elongated in time direction [49| . Let us assume that the geometry is 
indeed a prolate spheroid, for convenience we can imagine a two-dimensional ellipsoid or 
a surface of a cigar. For a given triangulation, we start the diffusion at the central slice, 
i.e. on the smallest equator. At the beginning, the diffusion front propagates in concentric 
circles, forming a connected shell, this is illustrated by the main top branch in Fig. 13.91 
There is only a small number of short outgrows. However, at some point the diffusion 
front reaches the antipode of the starting point and bifurcates into two disconnected fronts 
moving in opposite directions toward the tips of the spheroid. This point is denoted by 
the bifurcation radius Vbif- The bifurcation radius r^jj is visible on the graphs of diffusion 
trees in Fig. 13.91 as the point at which the main branch bifurcates into two branches. 
Here, r^f is identified with half of the shorter circumference lying in the central slice 
icv, and for A = 0.25,0.6 it is much shorter than the time extent. While decreasing 
A the bifurcated branches become less distinct, and eventually vanish, which suggests 
that the four-dimensional geometry approaches a spherical shape. Moreover, the Figures 
13.91 indicate that the four-dimensional spacetime geometry does not reveal a true fractal 
structure when looking at a shell decomposition of spacetime, as only short outgrows are 
present. It should be noted that the diffusion trees were obtained using a four-dimensional 
definition of the neighborhood. By contrast, in Chapter [5] we examine a similar shell 
decomposition within a spatial three-dimensional hypersurface and show that a fractal 
structure is indeed present. In the latter case we use a three-dimensional definition of the 
neighborhood of tetrahedra. 
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Figure 3.11: The asymmetry factor a between lengths of time-like and space-like edges 
of the triangulation as function of A for Kq = 2.2. The thin horizontal line corresponds 
to the boundary value 

Analyzing the volume distribution N^^i) (cf. Fig. 13. 8p and diffusion trees (cf. Fig. 13. 9p 
we argued that geometry of quantum spacetimes appearing in CDT resembles a perturbed 
four-dimensional elongated spheroid. To support this statement and the dependence on 
the coupling constant A, we use a more quantitative observable, namely (A^4(i,r)). In 
contrary to the diffusion trees, it does not refer to a single configuration, but is averaged 
over the ensemble of triangulations and initial points (for each configuration ca. 100 
starting simplices were chosen from the middle slab icv)- Fig- ETO] presents contour plots 
of the distribution {N4{i,r)) for A = 0.10,0.25 and 0.6. For large values of A > 0.15 it 
assumes a y-shape in the time-radius plane. This is an expected behavior for the elongated 
spheroid. We start the measurements of A^4(z,r) at r = which is a simplex in the initial 
slab icv which corresponds to a vertical line passing through the exact middle of the plot. 
Thus, the starting point N^^icvi^) matches to the bottom tips of the Figures r3. 101 When 
considering further shells, with increasing radius r, i.e. going upward in Figures 13.101 the 
shells spread over successive time slabs. However, at the bifurcation radius, the shells 
eventually reach the antipodal point completely covering the initial slab. In this case the 
bifurcation radius ri^j corresponds to the largest radius r for which {Ni{icv,f)) is larger 
than zero. For r > r^^f shells move away from the original slab, both r and i are growing 
causing the y-shape of the {N^^i, r)) plot and explaining the elongated shape of quantum 
geometries. As A decreases the V-shape becomes less pronounced (cf. Fig. I3.10p . with an 
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Figure 3.12: The average radial volume {Nj^ir)) for A = 0.1,0.25 and 0.6. {N4{r)) is 
defined as the number of simplices located exactly at discrete distance r, from an initial 
simplex lying in the middle slice icv average both over configurations and starting points. 



obvious interpretation that the shape of the spheroid becomes more and more spherical, 
with approximately equal extent in spatial and time directions. 

So far, the time extent of the Universe was described in terms of lattice spacing. 
However, in four-dimensional CDT the additional coupling constant A is related to the 
asymmetry factor a of the lengths of time-like and space-like links. It is defined by equation 
(|1.4|) . and affects the thickness of time slabs. The factor a is a function of Kq, A and the 
critical value of which can be obtained by inverting equations ()A.7p . When a grows, 
the slabs are getting wider. As seen in Fig. 13.111 this happens when we decrease A and 
approach the B-C phase transition line. Indeed, together with the decreasing time extent 
of geometries in terms of a number of slabs, the slabs themselves are getting wider and 
the two effects may compensate. 

The sphericity. Let us quantify what we mean by a sphericity of a triangulation. 
We introduce the sphericity factor s to measure how spherical is the averaged geometry, 

where the bifurcation radius rf,j/ corresponds to the radius r at which {N^^icvi^)) fe^s 
below some cut-off (here taken to be 4 to suppress the effect of large fluctuations), icv is 
the initial central slice. Such definition is in agreement with the average bifurcation radius 
present in diffusion trees. 
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Figure 3.13: Spliericity factor s for Kq = 2.2 and different values of coupling constant A. 

The distribution N4^(r), as defined at the beginning of this Section, denotes the number 
of four-simplices in a shell of radius r. Fig. 13.121 shows {N/^^r)) for various values of A. 
For the smallest values of A the peak is nicely symmetric and well approximated by 
A ■ sm^ (r/B), which corresponds to a sphere. For larger values of A, the elongation of the 
geometries is reflected in a large-r tail which cannot be fitted to A • sin^(r/S). 

The sum in the denominator of ()3.28p equals the total volume A'4 with the stalk 
cut away. The sphericity factor denotes the part of volume which is located up to the 
bifurcation radius. For a perfect sphere we have s = 1, as whole four- volume is located 
between the origin and antipodal point. For a very elongated object s is close to zero. 
Fig. 13.131 shows the sphericity factor s, for data averaged both over configurations and 
initial points, as a function of coupling constant A at fixed Kq = 2.2. When decreasing 
A and getting closer to phase B, the sphericity factor grows and approaches to s = 1, 
i.e. the shape of the average geometry changes from a prolate four-dimensional spheroid, 
elongated in time direction, and equal radii in spatial directions, into an almost spherical 
shape. This picture is consistent with what we observe for diffusion trees. 

To conclude, the dynamically emerged background geometry corresponds to the Eu- 
clidean de Sitter solution of the minisuperspace model. It may be very well described by 
a prolate four-dimensional spheroid for points deep inside phase C. When approaching 
phase B the Universe is getting narrower in terms of lattice spacings, finally reaching a 
spherical shape. However, this might be an illusion since one should take into account the 
change of slab width. Lastly, our more detailed investigation shows little or no evidence 
of a fractality when looking at a shell decomposition of the spacetime. 




Chapter 4 

Quantum fluctuations 



This chapter is based on the following publications: (i) J. Ambj0rn, A. Gorlich, J. Jurkie- 
wicz and R. Loll, "Planckian Birth of the Quantum de Sitter Universe", Phys. Rev. Lett. 
100, 091304 (2008); (ii) J. Ambj0rn, A. Gorlich, J. Jurkiewicz and R. Loll, "Nonpertur- 
bative quantum de Sitter universe", Phys. Rev. D 78 (2008) 063544; ('i'^'i) ^- Gorlich, 
"Background Geometry in 4D Causal Dynamical Triangulations" , Acta Phys. Pol. B 
39, 3343 (2008); (iv) J. Ambj0n, A. Gorlich, J. Jurkiewicz, R. Loll, J. Studnicki and 
T. Trzesniewski, "The Semiclassical limit of Causal Dynamical Triangulations" in prepa- 
ration. 



In the previous Chapter we gave the arguments that the dynamically emerging back- 
ground geometry can be obtained as a solution of the minisuperspace model. Here we 
will check if quantum fluctuations around the classical trajectory ()3.10p are also correctly 
described by the effective action (|3.9p . In this Chapter we investigate the properties of the 
semiclassical limit of the lattice approach. Nevertheless, it should be clearly stated that 
these considerations are truly non-perturbative, and take into account both a very impor- 
tant influence of the entropy factor, which does not depend on bare coupling constants, 
as well as the bare action (|1.5p . 

As before, we integrate out all degrees of freedom except the scale factor, and restrict 
our deliberations to three- volume, thereby reducing the problem to one-dimensional quan- 
tum mechanics. Based on numerical data obtained by computer simulations, we construct, 
within the semiclassical approximation, the effective action describing discrete spatial vol- 
ume N{i) and compare it with the minisuperspace action ()3.9p . The effective action comes 
into existence because of a subtle interplay between the entropy of configurations, which 
depends on the path integral measure, and the bare action (jl.Sp . Let us recall, that the key 
ingredient which allowed to perform a meaningful average of the volume over geometries, 
is the removal of the translational zero mode. The same is true in the case of quantum 
fluctuations of the spatial volume. 
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Let us denote the deviation of the three- volume N{i) from the expectation value N(i) 

by 

r]i = N{i) - N{i). 

Imitating the path integral approach to quantum mechanics, N{t) describes the position 
of a non-physical particle trajectory, giving a contribution to the partition function, where 
now t = i is the continuous time coordinate. Likewise, r](t) is a fluctuation from the 
classical trajectory N{t). In the semiclassical approximation, the spatial volume fluctu- 
ations r]{t) are described by a Hermitian Sturm-Liouville operator P{t), obtained by the 
quadratic expansion of the effective action around the classical trajectory 

S[N = N + r]]f^ S[N] + i y r]{t)P{t)rj{t)dt + 0{r]^). (4.1) 

In the CDT model, the time slicing carries with it the discreteness of time coordinate i, 
thus the Sturm-Liouville operator is substituted by a matrix P, and the expansion (j4.ip 
is 

5[iV = iV + 77] « S[N] + ^ E + (4-2) 

where the sum is performed over time slices i,j = 1 . . . T. 

The P matrix carries information about quantum fluctuations and may be extracted 
from numerical data obtained by computer simulations. In analogy to {N{i)) (cf. ()3.3p ). 
we measure the covariance matrix C of volume fluctuations using Monte Carlo techniques, 
described in detail in Chapter [6l 

k=l 

The covariance matrix C is also called the propagator. If the quadratic approximation de- 
scribes properly quantum fluctuations around the average N , the following direct relation 
between the propagator C and the matrix P is satisfied, 

= I / ??^r/,e-|Efc,..'?..P..'r,,, -Q^^^ ^ p-1^ 

s 

where Z denotes the normalization factor. Fig. 14.11 and Fig. 14.21 present the plot of the 
measured covariance matrix C and its inverse P for coupling constants Kq = 2.2, A = 0.6 
and the total volume Ntot = 80k. 

For the numerical convenience, measurements were performed only for triangulations 
with a fixed total volume Ntot = Yli=i^ (''')• This constraint imposes on the covariance 
matrix C the existence of a zero mode, preventing from obtaining the P matrix by a 
straightforward inversion of C. For each configuration the sum of volume fluctuations 
vanishes, rji = N{i) — N{i) = Ntot — Ntot = 0. Hence, the zero mode corresponds 
to a constant vector e*^. 
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In order to invert the matrix C we project it on a subspace orthogonal to the zero mode e 
and then perform the inversion. This can be achieved in the fohowing way. First, we add 
by hand to C the dyadic product of and e". This shifts the zero eigenvalue, which now 
equals one, and allows to invert the redefined matrix. Further the zero mode is recovered, 
by subtracting back the same dyadic product. The matrix P, defined as the inverse of 
C on a subspace orthogonal to the zero mode e^, is given by the following formula (cf. 
derivation of ()B.4|) ). 

P = (C + A)-i-A, A,, = eOeO = ^. 

The tentative procedure of dealing with the zero mode may be avoided by modifying 
the action by a quadratic term and allowing the system to perform measurements for an 
arbitrary total volume. Such correction removes the artificial zero mode and simplifies the 
analysis, but in fact does not affect results presented here. Together with the extended 
analysis taking into account also half-integer slices built of {3, 2} simplices, the new method 
will be addressed elsewhere [HT] . 

After calculating the covariance matrix C, we can get the empirical Sturm-Liouville 
operator P which can be compared with the predictions of the minisuperspace model. As 
may be seen from Fig. 14.21 the empirical P matrix has to a very good approximation a 
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Figure 4.2: The inverse propagator Pjj extracted from the numerical data for Kq = 2.2, 
A = 0.6 and Ntot = 80k. 



tridiagonal structure. Except for the diagonal and both sub-diagonals, the remaining part 
of the P matrix becomes indistinguishable from a numerical noise. Indeed, we will show 
that such structure is expected from the minisuperspace action (|3.9p . 

The tridiagonal form suggests that the effective action describing fluctuations of N{i) 
is quasi-local in time. 



Smz)] = /(iV(z), iV(z + 1)) + E U{N{i)). (4.4) 

i i 



The function f{x,y) describes the kinetic part of the effective action. It couples volumes 
of successive slices providing the non-zero subdiagonal elements of P. Because both C 
and P are symmetric matrices, the function / has to be symmetric in its arguments, 
f{x,y) = f{y,x). The potential part is described by function U{x) which contributes only 
to the diagonal. According to the expansion of the effective action (|4.2p the elements of 
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the Sturm-Liouville matrix P are given by 



N=N 



where 



dN{i)dN{j) 

/(i'i)(iV(i),iV(i-l)) J = i-l 

/(2.o)(iv(z),iV(i - 1)) +/(2.o)(iv(i),iV(i + 1)) + U"{N{i)) ,j = i 

/W)(iV(i),iV(i + l)) J = i + l 

, otherwise 



(4.5) 



Discrete minisuperspace action. In the next Sections, we determine functions 
f{x,y) and U{x), and show that the effective action (|4.4p corresponds to a discretization 
of the minisuperspace action ()3.9p up to an overall sign. Below we derive a discrete version 
of the minisuperspace action with reversed sign 



S[v] 



dr f a— + /3t;^/3 - 2Av 



(4.6) 



which later will be compared to the empirical action. We have incorporated the factor 
into constants a, (5 and A. The discretization procedure is not unique, but up to the 
order considered here, all discretizations are equivalent. We substitute the physical volume 
v{t) with the discrete volume A^(i), which as stated in Section [3.11 mav be treated as a 
continuous variable inside the hloh. The stalk region is governed by very strong lattice 
artifacts, and therefore is not reliably treated in the semiclassical approximation. The 
standard discretization of the time derivative is i) — )■ N{i + 1) — N{i), and the kinetic part 
is written as 

"■•2 (iv(f + l) -Ar(i))2 



a- 



91- 



N{i + 1)+N{i) 



while the potential part as 



Therefore, a discretized, dimensionless version of (j4.6p is given by 

(Ar(i + 1) _ iV(i))2 



S[N] = 



91- 



N{i + l)+N{i) 



+ 92N{i)^/^ -gsN{{), 



(4.7) 



which corresponds to 



x + y 



g3x. 



Further we show that the discrete effective action (|4.7p describes not only the average 
N{i) ()3.4p . what follows from the classical trajectory of equation ()4.6p . but indeed also 
the measured fluctuations r](i). 
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4.1 Decomposition of the Sturm-Liouville matrix 



From the formula (j4.4p follows that Sturm-Lfouville operator decomposes into the kinetic 
part P'^*", given by the function f{x, y), and the potential part PP°*, given by the function 
U{x), 



ykin 



_|_ ppot 



Only the kinetic part contributes to the sub-diagonal elements of the tridiagonal matrix 
P (cf. (|4.5p ). We make following assumptions about p'^*": 

• It should be a tridiagonal matrix. The square of the time derivative should couple with 
the preceding and following time steps. 

• It should be a symmetric matrix. Covariance matrix is symmetric, and so should be P. 

• The sum of elements in each row or column should be zero. A constant vector should 
be a zero mode eigenvector as it is invariant under time translations. 

Taking into account time-periodicity, this implies a following form of p'^*"^ 

/ kr + h -ki ... -kr \ 

-ki ki + k2 -k2 ... 

-k2 k2 + k3 -ks 



■tkin 



\ -kT 







(4.9) 



—kx-i kx^i + kr / 



The matrix Y* can be decomposed into parts linearly dependent on ki: 

T 
i=l 



(4.10) 



where X^*) is a matrix corresponding to the discretization of the second time derivative 



at a time t = i, 



In the case of the potential term, we have to get rid of the zero mode. The measured 
potential pP°* is a projection of the potential part of the full, unconstrained inverse prop- 
agator pP°* on a subspace orthogonal to e^. Assuming that the full potential matrix is 
diagonal 

PP°* = Diag({^za), 

and following equation ()B.3p derived in the Appendix B, we get the empirical propagator 
inverse 

^ 1 
PP°* = (1 - A) Diag({n,}) (1 - A) = ^n,Y», = e^eO = -, (4.11) 

i=l 

which can be decomposed into linear combination of matrices 

i - ^ 1 

T 
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Figure 4.3: Kinetic term: The directly measured expectation values N{i) (black line), 
compared to |^ (thick line) extracted from the measured covariance matrix C for Kq = 
2.2, A = 0.6 and various total volumes A'tot ranging from 20000 to 160000 simplices. The 
theoretical predication ^ = ^{N{i) + N{i + 1)) is realized with a very high accuracy. The 
value of gi is constant for all volumes A'tot- 



Here, 1 denotes the T x T unit matrix. The kinetic part remains unaffected by the 
projection, since AP'^*" = 0. 

The decomposition of the empirical matrix P into a kinetic and potential part, is done 
using the least square method. We find such values of {ki} and {ui}, for which the matrix 
pfcm _j_ ppot^ given by ()4.10p and (j4.1ip . is as close as possible to the empirical matrix P, 
i.e. we minimize the residual sum of squares 



RSS[{yfcJ, {ui}] =Tr P - (P'^^" + PP°* 



(4.12) 



We will omit details of the parameter fitting. In fact, we used the weighted least square 
method, giving higher weights to slices with larger volumes, i.e. lying inside the blob. 
Equation ()4.12p is quadratic in {ki} and {ui}, and the fitting boils down to calculating 
traces of products of matrices X^*) and Y'^-'^. 



60 



Quantum fluctuations 



0.25 



0.2 



0.15 



0.1 



0.05 




2000 4000 

m 



6000 



Figure 4.4: Tlie plot of 4 
Kq = 2.2, A = 0.6 and Ntot 



N{i+1) 



+ 4. 



N(i-1) 



2 as a function of N(i) for 



80k. If the minisuperspace action describes correctly 
fluctuations r](i), this function is proportional to the first derivative of the potential U'{x). 
The best fit a ■ x~'^ + h plotted with black line corresponds to c = 0.71 ± 0.04. The 
minisuperspace model predicts c = 1. 



4.2 Kinetic term 

In this Section we show that the fitted values of the kinetic term {A;^}, obtained by min- 
imizing residues (|4.12|) . are indeed in agreement with the kinetic part of the discrete 
minisuperspace action (|4.7p . For the minisuperspace model, equations (j4.5p and (|4.8p give 

8X7/ 



-Pii+i = -/(^'i)(x,y)=5i 



and in the zeroth order approximation, N{i) - 
the kinetic term 

{N{i) + N{i + l)f 



x = N{i), y = N{i + l), (4.13) 
N{i + 1), we have following behavior of 
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i(iV(.)+iV(i + l)). 



(4.14) 



8 N{i) ■ N{i + l) 

Fig. 14.31 presents the plot of gi/ki for the empirical values of ki and various total volumes 
A'tot- The theoretical fit (j4.14p agrees extremely well with the measured quantities. Addi- 
tionally, the effective coupling constant gi does not depend on Ntot in the margin of error. 
For Kq = 2.2, A = 0.6, we measured gi = 0.038 it 0.002. The kinetic part of quantum 
fluctuations is indeed described by the minisuperspace action ()4.6p . 
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This success in applying the minisuperspace model to Causal Dynamical Triangulations 
in four-dimensions, encourages us to check if we can say something about the potential 
term making use of the knowledge only about the kinetic part and the classical trajectory. 
The classical solution minimizes the action, 



dN{i) 



= 0. 

{Nii)=N{i)} 



Assuming that the minisuperspace model is valid, after inserting the function f{x, y) given 
by ()4.8p into the action ()4.4p we get an equation for the first derivative of the potential 
U{x) for arguments equal to the average spatial volumes, 

U'(N(i)) = 45i_^^^^±^l— - + 4gi ,-, - 251. 

Fig. 14.41 presents the evaluated derivative of the potential, for empirical values of average 
spatial volume N{i)- Again, the measurements agree within the limit of error with the 
prediction (|4.8p 

4.3 Potential term 

In this Section we directly show that values of the potential term {uj} extracted from the 
empirical inverse propagator P also agree with the minisuperspace model. Within this 
framework, following equations (j4.5p and (|4.8p . we expect that 

u, = U"{N{i)) = -^g2N{i)-'/\ (4.15) 

Fig. 14.51 shows the measured values of coefficients Ui extracted from the empirical 
matrix pp°*. Because of large statistical errors, it is not an easy task to determine Ui. As 
can be seen in Fig. 14.51 the plot is dominated by the stalk which is governed by lattice 
artifacts and can not be reliably treated semiclassically. The physically interesting region 
of large volumes corresponds to relatively small values of Ui as they are expected to falls 
as N{i)~^/^. Due to the projection on the space orthogonal to the zero mode, the blob 
region is affected by the huge contribution from the stalk. Moreover, analogically as in the 
ordinary path-integral approach to quantum mechanics when the time step approaches 
zero, in the continuum limit A'tot — ^ oo the potential term is sub-dominant w.r.t. the 
kinetic term for individual spacetime histories in the path integral. 

Nevertheless, due to sufficiently long Monte Carlo sample, the obtained results allowed 
us to confirm that indeed the formula (j4.15p is in agreement with the measurement. Fig. 
14.61 presents the measured coefficients — Uj as a function of the average three- volume N{i). 
The error bars shown on the plot, were estimated using Jackknife method. Such a form 
allows us to directly compare the potential coefficients with theoretical predictions —Ui oc 
N{i)~^^^. The selected range of N{i) corresponds to the bulk, and is highlighted in Fig. 
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Figure 4.5: Coefficients Ui of tlie potential part extracted from the empirical matrix P 
obtained in Monte Carlo measurements for Kq = 2.2, A = 0.6 and A^tot = 80k. The coef- 
ficients correspond to the second derivative of the potential Ui = U"{N{i)). The potential 
term is dominated by the stalk region. The physically interesting range, corresponding to 
the blob, is highlighted. 



14.51 The best fit of the form f{x) = a ■ x ^ to the empirical values function 
of N{i) gives c = —1.658 it 0.096. The measured exponent coefficient c is very close to 
the theoretical value c = —5/3. The fit f{x) = a ■ x~^^^, corresponding to the potential 
()4.8p is presented in Fig. 14.61 with a thin line. The agreement with the data is good, the 
potential part of the effective action is indeed given by U{x) = g2X^^^ — g^x. There is 
a very small residual constant term present in this fit, which presumably is due to the 
projection onto the space orthogonal to the zero mode. In view of the fact that its value is 
quite close to the noise level with our present statistics, we have simply chosen to ignore 
it in the remaining discussion. Apart from obtaining the correct power Uj oc 

7V-5/3(i), the 

coefficient in front of this term is also independent of A'tot- 

In summary, we conclude that the measurements of the covariance matrix of quantum 
volume fluctuations allowed us to reconstruct the discrete version of the minisuperspace 
action (|4.7p with a high precision. 



I I 




I I 



4.4 Flow of the gravitational constant 



63 



6e-06 



4e-06 



a 
I 



2e-06 



Fit av'^l^ 
Monte Carlo, 160k 




4000 



5000 



Figure 4.6: The extracted potential term Ui as a function of average volume N{i). The fit 
-—5/3 

C2V^ presents the behavior expected for the minisuperspace model. The visible points 
correspond to the blob region. 

4.4 Flow of the gravitational constant 

In the previous Section we gave a direct numerical evidence that quantum fluctuations of 
three-volume are very accurately described by the discrete, dimensionless effective action 



S[N{i)] 



1^91 



2N{i) 



+ 92 



(4.16) 



where we omitted the cosmological constant term since during the measurements the total 
volume A'tot was fixed. This action comes out as a discretization of the minisuperspace 
action (|3.9p with the opposite sign which solves the problem of unboundedness. Let us 
note, that it is justified to use the semiclassical approximation as the distribution of spatial 
volumes N{i) in the bulk (cf. Fig. 13. 2p is given by Gaussian fluctuations around the mean. 

In this Section we establish the relationship between the discrete action (j4.16p and its 
continuum counterpart (|3.9p , which will allow us to connect the effective coupling constant 
gi with the gravitational constant coupling constant G. At first, let us rewrite the above 



action in terms of the normalized volume n{t). Using equations N{i) 



3/4 



and At = N^^^^ (fXTB . where t 



S[n{t)] = 



n{t) ([3^2]) 



At • i, the discrete action can be written as 
1 / n{t + At) - n{t) 



91 



2n{t) 



At 



+ 92 n^'\t) 
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and following (j3.13p . we can represent it in a form of an integral over the continuum time 
t 



S[n{t)] = / dt 



1 i/s/ N 



(4.17) 



As stated in the previous Sections the effective coupling constant gi does not depend on 
A'tot when the bare coupling constants are fixed, the same is true for the classical trajectory 
n{t). Hence, the amplitude of fluctuations of n(t) scales with the total discrete volume A'tot 
as ((5n(t))^) oc ^^"^N^^J:^^ and they vanish in the continuum limit A'tot oo, while the 
background n{t) stays fixed. Regarding Fig. 13. 4| it means that the height of the curve will 
grow as N^^t^ superimposed fiuctuations will only grow as , and for fixed bare 

coupling constants the relative size of fluctuations will go to zero in the infinite-volume 
limit. Therefore, in order to obtain a proper continuum limit with a finite physical four- 
volume V4 (X Nx,ot<i^ = const and finite fiuctuations \/T{^v{Typ^ oc g^^ ^^"^ N^l^a^ = const 
(cf. (j3.20p ). one also has to properly tune the hare coupling constants so that the effective 
coupling constant satisfies gia~'^ = const while taking A'tot 00 and a — t- 0. 

Further, we rewrite the action (|4.17p in terms of the physical volume v{t). Using 
formulas ()3.19p and ()3.20p . we recover a continuum action back from the discrete action 



Slv(r)] = ^ I dr 



v{r) 



+ 92 v^'\t) 



(4.18) 



This is the second most important result obtained in this dissertation. Not only the 
average spatial volume v{t) obtained in the Causal Dynamical Triangulations frameworks 
is given by the classical trajectory of the minisuperspace model, but also the fluctuations 
of v[t) around the background geometry v{t) are described by the minisuperspace model 
up to the overall sign. 

It is natural to identify the coupling constant G multiplying the effective action for 
the scale factor (|3.9p with the Newton's gravitational constant G. The effective action 
describing our computer-generated data (|4.18p can be directly compared with the min- 
isuperspace action p.9p . and using equations (j3.2ip and ()3.22p . we get following relations 
between the gravitational constant G and the effective constant gi [351 HB], 



24^91 fftt gi 3^/6 

This result is consistent with our earlier conclusion, in order to keep the physical constant 
G fixed, when taking the continuum limit a — t- one has to tune the effective coupling 
constant gi oc a^. This means that in terms of the lattice volume N{i) fiuctuations should 
diverge, and this happens when we approach a second or higher-order transition line. 
Therefore it is important to determine the order of transition. If the B-C transition would 
be of the first-order we may hope that the triple point is of higher order, which is often true, 
and obtain the proper continuum limit while approaching this point. Fig. 14.71 shows the 
measured effective coupling constant gi for various values of A. Indeed when we approach 
the B-C transition line gi tends to zero. However, it is not so easy to find the exact path 
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Figure 4.7: The measured effective coupHng constant gi as a function of bare coupling 
constant A for Kq = 2.2. The B-C transition point is located c.a. at A^"* = 0.05. 
When approaching phase B from phase C the coupling constant gi diminishes and the 
fluctuations grow, as expected when reaching phase transition point. 

in Kq — A coupling constant plane, as formula (j4.19p depends also on other parameters, 
namely the effective four- volume shared by one spatial tetrahedron 6*4(0;) and the width 
B of the distribution n(t), and it will require further extensive computer-experiments. 

The Universe size. Using relation ()4.19p we can express the cut-off length a in terms 
of the Planck length, and thus estimate the size of the Universe generated in computer 
simulations. Let us recall that in natural units G = £pj. For the bare coupling constants 
Ko = 2.2, A = 0.6 we measured the quantities: Kf''^ = 0.922,^ = ^ = 1.30, a = 
0.5858, C4 = 0.0317,(^1 = 0.038, which results in a ~ 1.9ipi and the linear size irR of 
the universe built from 160000 simplices is about 20^p/. The quantum de Sitter universes 
studied here are therefore quite small, and quantum fluctuations around their average 
shape are large (cf. (|3.4|) ). Surprisingly, the semiclassical minisuperspace formulation 
gives an adequate description of the measured data, at least for the volume profile. 

To conclude, based on the numerical evidence we reconstructed the discrete effective 
action describing quantum fluctuations of the three- volume N{i). This allowed us to 
identify the effective action with the discretization of the minisuperspace action with a 
negative sign. We also derived a relation between the physical gravitational constant and 
the effective coupling constants appearing in the model, giving a recipe of how to obtain 
a meaningful continuum limit. 



Chapter 5 
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This chapter is based on the publication: J. Ambj0rn, A. Gorlich, J. Jurkiewicz and 
R. Loll, "Geometry of the quantum universe", Phys. Lett. B 690, 420 (2010). 

After investigating in Chapter [3] the four-dimensional geometry of the quantum Uni- 
verse, the next step is to look deeper into the geometry of spatial slices. A spatial slice is 
a leaf of the imposed global proper-time foliation and is labeled by a discrete time index 
i. Each slice is built of equilateral spatial tetrahedra. A hypersurface of this kind is a 
three-dimensional triangulation, more precisely, a piecewise linear manifold of topology 
. However, it does not mean that the geometry of slices is close to the geometry of a 
three-dimensional sphere. 

5.1 Hausdorff dimension 

Since in this Chapter we consider each slice as a separate object, let us denote the number 
of tetrahedra building a slice i by the discrete three- volume = N{i). A basic observable 
defined on a slice, is the number of tetrahedra n(r, zq) at a three-dimensional distance r 
from some initial tetrahedron ig- A three-dimensional distance is based on the definition 
of tetrahedra neighborhood within one slice, while a four-dimensional distance is based 
on the definition of simplices neighborhood (cf. Section [3. 4p . This difference introduces a 
significant change in the shell decomposition picture. At distance r = only the initial 
tetrahedron is counted and n(0,io) = 1- Next, we move out by one step visiting all 
four neighbors of io, thus n{l, zq) = 4. Moving out by a further step, there will be n(2, zq) 
tetrahedra at a distance r = 2, and so on until all tetrahedra are visited, with the restriction 
that one tetrahedron is visited only once. For such definition, n(r, io) corresponds to an 
area of the shell of radius r. Summing up the area over all shells gives the discrete volume 
of a slice ns. 
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Figure 5.1: The cloud visible on the plot presents the set of measured points (713, f). The 
darker pixel, the more events were observed in given range. To clarify the picture, green 
points correspond to an average of 100 successive measurements sorted by n^- The blue 
line represents the fitted curve f oc 71-3^"''' using linear regression, = 2.984 it 0.024. The 
plot is in log-log scale. 



Further, n(r) denotes the average of n(r, zq) over all 723 initial tetrahedra iq, 

n(r) = — y~] n{r,io). 

10 = 1 

In this Section we investigate scaling properties of n(r) with respect to the slice volume 

The Hausdorff dimension of spatial slices may be measured, using the scaling proper- 
ties, in a following way [30]. For a given slice of volume n^, we define the average linear 
extent as 

f = — > r • n(r). 
^3 V 

For each slice we obtain a data pair {723, f}. Such pairs are then gathered for all slices of 
a number of Monte Carlo (MC) configurations. A plot off all points belonging to some 
MC sample, shown in Fig. 15.11 for = 160A;, looks like a cloud. To make the average 
behavior more visible, the number of data points is reduced. The results are sorted by 
7i3 and averaged over a sequence of 100 consecutive points (in contrary to [30] where the 
sample is sorted by f, because fluctuations of f are larger than n^, the approach used 
here gives smaller finite size effects and better linear behavior of the log-log plot.) For 
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Figure 5.2: Scaled values of the radius ''r and shell area '^n{r) for = 3. 

Data points for various values of slice volume ns overlap. The gray strip plots the scaled 
radial volume averaged over all data points. The fitted curve ar^ (black line) presents the 
small r behavior and is consistent with Hausdorff dimension dh = 3. Measurements were 
performed at Kq = 2.2 and A = 0.6. 



Hausdorff dimension d/^, the expected relation between the slice volume 713 and the average 
linear extent r is following 

n3 0cf"'\ (5.1) 

The reduced data points (green dots) together with the relation (|5.1|) fitted to the full data 
are presented in Fig. 15.11 The fitted value of the Hausdorff dimension is df^ = 2.984ib0.024. 

Similar results can be obtained by a direct comparison of the scaled with volume nj, 
values of the radial volume n(r). First, for a large number of Monte Carlo configurations, 
slices with the same volume ns (more or less few tetrahedra) are collected into groups. 
The average radial volume n(r) within a group n-^ is denoted as {n{r))nz- Next, the radius 
r and the average volume {n{r))ns, are scaled and normalized in the following way 

(r, (n(r))) ^ (n~"'''^r, n'^+^Z'^'^ (n(r))) . (5.2) 

Analogically as in Section [3.41 where we calculated the scaling dimension du for N{i) (cf. 
Fig. 13.51 and 13. 6p . we define the error of the overlap of the scaled points and find such 
value of dh which minimizes the dispersion. The best fit is obtained for dh = 2.94 it 0.05. 
Fig. 15.21 presents the measured values of {n{r)) scaled according to ()5.2p with dh = 3 and 
for various values of na between 1000 and 4000 tetrahedra. 
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The HausdorfF dimension may be also estimated from the smah r behavior of 
(f^('"))n3- The value dh determines the dependence of the volume of a ball on its radius r, 
Vol[-B(r)] oc r'^'^ or equivalently for an area of a sphere 

Area[S'(r)] oc r'^'^^^, when r — t- 0. 

Because n(r) is interpreted as the area of a sphere with a radius r, the expected behavior 
of (n(r))„3 for small r is (n(r))„3 oc r'^^-^^ Assuming that {n[r))n^ oc r° for small r, 
a = dfi — I is the only exponent consistent with the scaling (|5.2p . Fig. 15.21 shows that 
dh = 3 is in agreement with the measured data. This method of estimating the Hausdorff 
dimension, is however more sensitive on the fit region than studying the scaling properties. 

To conclude, all methods of determining the Hausdorff dimension of spatial slices dh 
give consistent results. The measured value is with a high accuracy equal to dh = 3. 
Moreover, the measured value of dh is independent of the coupling constants Kq and A, 
as long as we stay well inside the phase C. This results is true if we consider the ensemble 
average of the slice geometry. As we shall see in the next Sections, this result does not 
mean that individual spatial slices resemble a smooth three-dimensional geometry. 



5.2 Spectral dimension 

In this Section we determine spectral dimension ds by investigating diffusion process on 
triangulated slices dominating in the CDT partition function in phase C. We measure the 
spectral dimension in the same way as described in Section [3.41 The probability of finding 
a diffusing particle in tetrahedron i after a diffusion time a and starting at tetrahedron io 
is given by the probability density p{i,iQ;a). The discrete diffusion equation, describing 
the evolution of the probability density, is an analogue of the equation ()3.26p . namely 

p(z,zo;cr + 1) = ^^p(j,io;cr), (5.3) 

where the sum is over all tetrahedra j adjacent to i. For a starting tetrahedron io, chosen at 
random, we set the initial condition p{i, iq; o" = 0) = (Jjjp. By iterating the diffusion equa- 
tion, we calculate the return probability P{cr, zq) = p{io, io', cr) for successive discrete diffu- 
sion steps a. Further, we compute the average return probability P{(t) = ((-P(<t, io))io)MC 
by averaging over initial points and configurations. For each configuration we consider 
only the central slice. Spectral dimension dg is obtained from the return probability using 
the definition ^^^2^ 

. ,5.4) 

dlogfj 

In the case of Euclidean space M*^, the spectral dimension and Hausdorff dimension are 
equal to the topological dimension dg = dh = d. For a three-sphere S^, spectral dimension 
dg is equal 3 for short diffusion times. Because of the finite volume, for longer times the 
zero mode of the Laplacian will dominate and dg will tend to zero. Fig. 15.31 shows values 
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Figure 5.3: Spectral dimension dg of spatial slices as a function of diffusion time a. For 
short diffusion times, a split for even (blue) and odd (teal) values of a is observed arising 
from the discrete structure. The measured values of dg converge to thin line corresponding 
to ds = 1.5. 

of the spectral dimension dg as a function of the diffusion time a determined by numerical 
simulation using the definition ()5.4p for a randomly chosen typical configurations in phase 
C. Due to the discrete lattice structure, for small values of a a split for even and odd 
diffusion times is observed. Because of finite volumes of slices, for very large a, dg falls 
down to zero. For the intermediate region, there is a plateau of the spectral dimension at 
ds w 1.5. 

5.3 The fractal structure of spatial slices 

The measured value of the spectral dimension of spatial slices dg ~ 1.5 is significantly 
smaller than the Hausdorff dimension d/j ~ 3. The difference between dh and dg is an 
indication of a fractal nature of slices. We will now identify this fractal structure in a 
more direct way. Three-dimensional spatial slices demonstrate a very similar structure to 
triangulated universes present in three-dimensional Euclidean Dynamical Triangulations 
(EDT). It is known, that in three-dimensional EDT so-called baby Universes separated by 
minimal necks are observed |641 165] . Advocated by the presence of such objects in the 
Euclidean theory, we look for their existence inside spatial slices. A minimal neck de- 
notes a set of exactly four triangles forming a tetrahedron (i.e. each triangle is connected 
to each other), which is not present in the triangulation. It corresponds to the smallest 
nontrivial boundary of a three-dimensional simplicial manifold built of tetrahedra. The 
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Figure 5.4: A visualization of a minimal neck on a two-dimensional triangulation. The 
minimal neck is marked with a thick black line consisting of three segments creating a 
triangle which itself is not a part of the triangulation. Cutting the triangulation along the 
minimal neck will separate it into two components: blue and red. 

minimal necks provide the three-dimensional triangulation with a tree-structure. A sim- 
plicial manifold cut along a minimal neck, by removing connections between tetrahedra 
which share the triangular faces belonging to the minimal neck, is decomposed into two 
disconnected parts. The two components form a simplicial manifold with a boundary and 
both of them contain more than one tetrahedron (in fact at least four). By closing the 
boundaries with tetrahedra we make the manifolds into two triangulation of a topology 
S^. Fig. 15.41 illustrates a minimal neck in the case of a two-dimensional triangulation. In 
two-dimensions a minimal neck is created of three links forming a missing triangle and is 
drawn with a thick black line, it separates the triangulation into a blue and red part. 

It is natural to visualize this structure in the form of a tree graph. The procedure of 
creating such tree is following: 

• We find all minimal necks, i.e. sets of four triangles forming a tetrahedron which 
does not belong to the triangulation. 

• Each triangle is a common face of two tetrahedra. Ungluing the tetrahedra sharing 
elements of a minimal neck decomposes the triangulation into two parts. 

• We cut the triangulation along all minimal necks, separating the slice into many 
disconnected pieces. We represent each piece as a vertex. To each vertex, we can 
assign the number of tetrahedra building a corresponding piece. This number is 
called a vertex volume. 

• We add an edge linking two vertices if they share a common minimal neck. Because 
the slice has a topology of a three-sphere, by definition the set of vertices and edges 
forms a connected tree graph, called a minhu tree. 

Fig. 15.51 presents two minhu trees created by the above procedure and corresponding to 
two randomly chosen successive slices. Considering one graph, we see that it is indeed 
a tree, but what is most important it has a fractal structure (regarding for finiteness). 
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Figure 5.5: The tree structures of minimal necks, siiown liere for two consecutive spatial 
hypersurfaces, reveal a fractal nature. The left tree corresponds to a demonstration slice 
i, while the right to slice i + 1. Some branches in tree i are marked with a chosen color, 
which also taints close vertices in adjacent tree. The colors identify neighbor vertices in 
two hypersurfaces. The diameter of vertices grow with their volume. 

Very often, but at random places, a branch bifurcates into two or more branches. Most 
probably, when the size of the slice grows to infinity, one would observe a real fractal 
structure of an infinite branched polymer. 

The definition of a minbu tree gives very similar results as the definition of diff'usion 
trees used in case of spatial slices (cf. Section [3^ . The full four-dimensional geometries 
did not exhibit any fractality, four-dimensional analogues of minimal necks appeared very 
seldom. In contrast, if we had created a diff'usion tree for a spatial slice, we would observe 
the same fractal structure as presented in Fig. 15. 5| with even more bifurcation points. The 
advantage of a minbu tree is that it does not introduce any specified initial tetrahedra. 
The graph in Fig. 15.51 is centered on the edge which divides the tree into two most equal 
parts with respect to the three-volume. 

The fractal structure suggest that leaves of the imposed time foliation have no semi- 
classical meaning. One has to perform an average of observables over all triangulations 
as we did in previous Sections. On the other hand, we checked that three-dimensional 
distance hierarchy reflected by the minbu tree (every path connecting two tetrahedra in 
different vertices must cross the appropriate minimal necks) , is not an artifact of the im- 
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posed foliation and persist in the full, four-dimensional embedding. There are no drastic 
shortcuts through the four-dimensional simplicial complex. 

This is corroborated by studying if the fractal structure survives between one time slice 
and the next. Fig. 15.51 illustrates minbu trees for two consecutive slices i and i + Using 
the full four-dimensional triangulation, we introduce the notion of neighbors on vertices of 
two trees. A vertex of a tree z is a neighbor of the vertex of a tree i + if there exists a pair 
of tetrahedra respectively shared by each vertex between which four-dimensional distance 
is 4, i.e. the shortest path between two slices passing through simplices {4,1}, {3,2}, 
{2, 3} and finally {1, 4}. Intuitively this means that two vertices are as close as possible in 
the four-dimensional simplicial manifold. One vertex may have many neighbors and the 
relation is symmetric. This notion allows us to study if the tree structure propagates over 
successive slices, i.e. if it survives from one slice to the next or if there is no correlation. 
For purposes of illustration, as shown in Fig. 15.51 we have marked all vertices of a given 
tree branch at time i in a chosen color. For each vertex we then find its neighboring 
vertices in time slice i -|- 1, and mark them with the same color, allowing also for mixing 
of colors. Clearly one can explicitly identify branches in the two slices, which proves that 
the fractal structure is correlated and not totally independent of each other. Nevertheless, 
the branches may significantly change their shape. The relation with further slices decays 
fast as the time step separation grows. So far, we have only this qualitative notion at our 
disposal. It is a challenge to define more quantitative methods measuring the correlation 
of successive trees, i.e. an analog of scalar product acting in the space of trees based on a 
four-dimensional definition of distance. 

We conclude that the measurements of Hausdorff dimension evidently show that = 
3. The spectral dimension ds is equal only half of this value, ds = 1.5. Such behavior 
indicates that the geometry of spatial slices is completely different from a smooth geometry 
of a three-sphere. Indeed, direct measurements on individual slices prove that they are 
equipped with a fractal structure similar to branched polymers. Quantum fluctuations can 
not be described by small Gaussian fluctuations around some average geometry, suggesting 
that the specifled hypersurfaces of constant time are not physical. Still, the tree structure 
propagates over consecutive slices. 
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The main tool used to obtain results within the framework of four- dimensional Causal 
Dynamical Triangulations are numerical Monte Carlo simulations. In CDT the partition 
function ()1.2p is regularized by a discrete sum ()1.9p over causal triangulations T. Such 
approach allows to generate Monte Carlo triangulations and calculate expectation values 
of observables according to (|1.10p . Therefore, we need to handle elements of the set of 
labeled causal triangulations T. In this Chapter we introduce software solutions and data 
structures used to store information about objects describing a triangulation. Later, we 
describe details of the implementation of the Monte Carlo algorithm, especially the Monte 
Carlo moves. 

Experience has shown that simulations are limited mainly by the available CPU time, 
rather than by the amount of the used memory. Therefore, the implementation of the 
algorithm uses a maximal information philosophy. The program stores in memory a full 
information about points, links, triangles and simplices. This approach allows for quicker 
access to the needed information at the expense of the increased memory demand. 

All results presented in this work were obtained using a computer program written in 
C. Therefore, we use the C language syntax to describe data structures. 

6.1 Parametrization of the manifold 

To perform the average (jl.lOp the Monte Carlo algorithm has to generate configurations 
according to the probability distribution on an ensemble of causal triangulations T. We 
have to define the way we code a simplicial complex and store it in the computer memory. 
To increase the acceptance of Monte Carlo moves, and thus the efficiency of the algorithm, 
we introduce many lists, enumerating the Lorentzian structures. 

Once we have labeled vertices of a triangulation and set its topology, the simplicial 
complex is given by the adjacency relations between four-simplices. Four-simplices, which 
are also labeled, share five vertices. The connections between simplices induce the adja- 
cency relations on all sub-simplices to which we also attach labels. Each label is an integer 
number between 1 and some maximal value. A label uniquely identifies sub-simplices ap- 
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pearing in the triangulation, and instead of writing "label of an object" we shall use a 
shortcut "object". The topology x is ensured by construction. As we shall see, we 
start with a minimal configuration with the designated topology and introduce only local 
Monte Carlo moves which do not spoil the topology. Therefore, to perform simulations 
the program has to code such objects as points, edges, triangles and four-simplices. It 
is not necessary to store information about tetrahedra (3-simplices) which if needed may 
be easily extracted. In this Section we describe most of the data structures used in this 
context. 

Points 

The points are labeled by a dynamical sequence of labels. The number of all accessible 
point labels is denoted as nnO. Therefore possible vertex labels are 1 ... nnO. The 
number of used point labels numO is equal to the number of vertices in the triangulation Nq. 
Some Monte Carlo moves are possible to be realized only on a local structure containing 
a point with a coordination number equal 8. Coordination number of a sub-simplex is 
denoted as the number of four-simplices sharing a given sub-simplex. The list of used 
labels is stored in the array listO which has a following structure: 

• Positions 1 ... numO are the used labels: 

— Positions 1 ... numSO are the labels of points with coordination 8 (used in 
move 4). 

— Positions numSO + 1 . . . numO are the rest of used (labels of) points. 

• Positions numO + 1 . . . nnO are the unused labels. 

The properties of all points are stored in the array p indexed by integer point labels: 
point p [nnO] ; 

The structure point (declared in the file types. h) holds information about a vertex: 

typedef struct 
{ 



int 


i; 


/* 


Index in listO 


*/ 


int 


n; 


/* 


Coordination number 


*/ 


int 


1; 


/* 


a random link which shares this vertex 


*/ 


int 


t; 


/* 


a random triangle which shares this vertex 


*/ 


int 


s; 


/* 


a random simplex which shares this vertex 


*/ 


int 


time ; 


/* 


time coordinate, number of the time slice 


*/ 


} point; 











The fields of the structure point have a following interpretation: 

• Field i contains index of this vertex in the array listO, 
i.e. listO [p [num] . i] = num. 
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• The coordination number n denotes the number of simphces, which share this vertex, 
n is not less than 8. 

• Field 1 is a label of a (random, among many) link sharing this vertex. 

• Field t is a label of a (random, among many) triangle sharing this vertex. 

• Field s is a label of a (random, among many) simplex that shares this vertex (used 
in move 4). 

• Field time stores the integer time coordinate of a point, i.e. the number of a slice to 
which the point belongs. Two points may be connected by a link only if their time 
difference is not larger than one. 

Links 

Links are labeled by a dynamical sequence of labels. The list of used labels is in the array 
listl with a following structure, analogical to array listO: 

• Positions 1 ... numl are the used labels: 



— Positions 1 ... num61s are spatial links with a coordination 6 (used in move 



— Positions num6 Is + 1 ... num61s + iium41t (next num41t) are time links with 
coordination 4 (used in move 2). 

— Positions num61s + num41t + 1 ... numl are the rest of used links. 
• Positions numl + 1 . . . nnl are the unused labels. 

The data about all links are stored in array 1 indexed by links labels: 
digon 1 [nnl] ; 

The structure digon (defined in types. h) holds information about a link: 



5). 



typedef struct 
{ 

int p [2] ; 



/* 
/* 
/* 
/* 
/* 
/* 



Labels of endpoints */ 

Coordination number */ 

Simplex which shares this link */ 

Position in listl */ 

Hash code of this link */ 



int n; 
int s ; 
int i ; 
int h; 
int n2; 



Next link in a hash table list 



*/ 



> digon; 



The meaning of the fields is as follows: 
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• The labels p of the endpoints. The labels are sorted: p[0] < p[l] . 

• The coordination number n denotes the number of simplices, that share this link, n 
is at least 4 for time links, and at least 6 for spatial links. 

• Field s is a label of one (random, among many) four-simplex containing a link with 
given label. 

• Field i contains the position of the link in the array listl: 
listl [1 [num] . i] = num. 

• Field h contains the hash code for the link: 

1 [num] . h = hlc (1 [num] . p [0] , 1 [num] . p [0] ) . 

• Field n2 contains the label of the next link in the singly-linked list of links with the 
same hash-code h. 

Hash map 

In move 2 a time link with a coordination 4 is created, while in move 5 a spatial link with 
coordination 5 is created. It might happen that a link with the same endpoints already 
exists. In this case we have to abort the move. 

In order to perform a fast search of a link given the labels of its endpoints, the hash 
table is used. Knowing the labels of vertices p[0] < p[l] we calculate the hash code h 
using the hash function hlc 

h = hlc(p[0], p[l]) = ((p[0] « 10) xor p[l]) 

defined in file core/hash. h. The first element of a singly-linked list of links with the same 
hash code h is given by hit [h] . The next element in the list is accessed via the field n2 
in the structure digon. To find if a given link exists, we have to look through this list - 
very short on average. 

Triangles 

The triangles are labeled by a dynamical sequence of labels stored in the array list 2: 

• Positions 1 ... num32t are time triangles with a coordination 3 (used in move 3). 
There are no time triangles with coordination less than 3. 

• Positions num32t + 1 ... num32t + num42s are spatial triangles with a coordina- 
tion 4 (used in move 5). There are no spatial triangles with a coordination less than 



4. 



Positions num32t + num42s + 1 



num2 are the rest of used links. 



Positions num2 + 1 



nn2 are the unused labels. 



6.1 Parametrization of the manifold 
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The properties of all triangles are stored in the array t indexed by a triangle label: 
triangle t [nn2] ; 

The structure triangle (declared in types. h) holds information about a triangle: 

typedef struct 
{ 

int p[3]; /* Labels of vertices */ 
int n; /* Coordination number */ 

int s; /* simplex which shares this triangle */ 

int i; /* Position in list2 */ 

int h; /* Hash code of this triangle */ 

int n2; /* Next triangle in hash table list */ 
} triangle; 

• Field p stores sorted labels of vertices: p[0] < p[l] < p [2] . 

• Field n is a coordination number of a triangle, i.e. athe number of simplices, that 
share this triangle. It is not less than 3. 

• Field s is a label of a simplex which shares this triangle (used in move 3). 

• Field i contains index of this triangle in the array list 2. 

• Field h contains the hash code of a triangle: 

t[num].h = htc (t [num] .p [0] , t[num].p[0], t [num] .p [2] ) . 

• Field n2 contains the label of the next triangle in the singly-linked list of triangles 
with the same hash-code h. 



Hash map 

In move 3 a time triangle with a coordination 3 is created, while in move 5 a spatial 
triangle with a coordination 4 is created. It might happen that this triangle already exists 
(a triangle with the same vertices). In this case we have to abort the move. 

To perform a fast check if a triangle with given vertices exists, a hash table is used. 
Given the labels of vertices p[0] < p[l] < p[2] we calculate the hash code h using the 
hash function htc 

h = htc(p[0], p[l], p[2]) = (p[0] * 256) xor (p[l] * 16) xor p[2] 

defined in a file core/hash. h. The first element of a singly-linked list of triangles with 
the same hash code h is given by htt [h] . The next element in the list is accessed via the 
field n2 in the structure triangle. Similarly as for links, to find if a given triangle exists, 
we have to look through this list - very short on average. 
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Simplices 

In CDT, there are two types of simplices: 

• type {4, 1} and {1, 4} witli four vertices in one time layer, 

• type {3, 2} and {2, 3} with three vertices in one time layer. 

The list of simplex labels is stored in the array list4 with a following order: 

• Positions 1 ... num444 are simplices of type {4, 1} and {1,4}. 

• Positions num444 + 1 ... num4 are simplices of type {3,2} and {2,3}. 

• Positions num4 +1 ... nn4 are the unused labels. 

The properties of all simplices are stored in the array s indexed by a simplex label: 
pentachoron s [nn4] ; 

The structure pentachoron (declared in types. h) holds information about a simplex: 

typedef struct 
{ 



int 


p[5] ; 


/* 


Labels of 


vertices */ 


int 


s[5] ; 


/* 


Labels of 


adjacent 4-simplices */ 


int 


1[10] ; 


/* 


Labels of 


links shared by simplex */ 


int 


t[10] ; 


/* 


Labels of 


triangles shared by simplex 


int 


i; 


/* 


Position : 


in list4 */ 


int 


rez [5] ; 


/* 


Reserved. 


Aligned to 32 */ 



} pentachoron; 

The fields of the structure pentachoron store a following information: 

• Field p stores labels of vertices of a simplex. 

• Field s stores labels of the neighboring simplices. Simplex s [i] is the neighbor 
opposite to p [i] . 

• Field 1 stores labels of links shared by a simplex. The links in the list 1 are sorted 
w.r.t. the first endpoint, and then w.r.t. the second one, where the ordering of the 
endpoints is given by p [i] . 

• Field t stores labels of triangles shared by a simplex. The triangles sorted w.r.t. the 
first, second than third endpoint where the ordering of endpoints is given by p [i] . 

• Field i is the index of a simplex in the array list4. 

On simplices, we finish the description of the data structures used to store full information 
about the triangulation. The general idea of the Monte Carlo program boils down to 
initialize above structures and sequential calling procedures that modify data in accordance 
with Monte Carlo moves. 



6.2 Monte Carlo Simulations 
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6.2 Monte Carlo Simulations 

The idea behind the Monte Carlo simulations is to approximate the sum appearing in the 
partition function ()1.9p by a sum over a finite number of Monte Carlo configurations. More 
precisely, given an observable 0\T], e.g. a distribution of a three- volume N{i) of spatial 
slices, we would like to calculate its expectation value {0\T\). In the CDT framework the 
geometries [g\ are already restricted to a discrete set of simplicial manifolds T and 

ref ^' 



where the partition function Z is given by eq. ()1.9p with the action (jl.Sp and the T 
corresponds to the labeling of T. 

Monte Carlo simulations generate a finite set of configurations {T^'^\ ■ ■ ■ ,T^^^} and 
allow to approximate the average (|6.ip by a summation over it, 

1 ^ 

^^[^)^:^E^[^^^^]- (6.2) 

i=l 

Let us notice, that no factor -T-^:^,e~^^ is needed, since configurations are generated 

No \T\ ! 

according to the probability distribution -P[T] = ^ |^, e"'^["^, which means that more 

probable geometries will more likely appear in the set {T^^\ ■ ■ ■ , T^^^}- We implement the 
Metropolis-Hastings algorithm to generate a Markov Chain using Monte Carlo methods 
[671 [68]. We define and implement local transformations, called Monte Carlo moves, which 
modify a configuration A. The Monte Carlo move applied to a triangulation A transforms 
it into a triangulation B. A series of such moves defines a Markov chain, the next configu- 
ration depends only on the last and on the applied move. After one step configurations A 
and B are almost identical, but after sufficiently large number of moves separating them, 
the initial and final configurations are independent. If the moves are chosen randomly, the 
algorithm performs a random walk in the phase-space of configurations, i.e. causal trian- 
gulations T. The moves are however not chosen completely at random. The moves and 
their weights, i.e. probability of being realized, must fulfill some constraints in order to 
properly probe the ensemble of causal triangulation T according to ()6.ip . The algorithm, 
after satisfying the conditions described below, ensures that it samples the configuration 
space according to a probability distribution 

The Monte Carlo moves must fulfill following conditions: 

Ergodicity. The moves generating the Markov chain have to be ergodic. It means that it 
is possible to reach any topologically equivalent causal triangulation by a sequence of such 
moves. Our scheme of assembling an appropriate set of elementary moves is to first select 
moves that are ergodic within the spatial slices. This is clearly a necessary condition for 
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ergodicity in the full triangulations. Namely, those are moves 4 and 5 together with their 
inverse moves 0. These moves act on spatial slices like three-dimensional Pachner moves 
which are equivalent to the primal Alexander moves and thus proved to be ergodic and 
preserving the 5^ topology of the space [70 1 169} ITT]. 

Then we supplement them by moves 2 and 3 acting within two spatial slices. They do 
not affect connections between tetrahedra building the spatial slices and can be understood 
as a Lorentzian variant, i.e. combination, of the four-dimensional Pachner moves, in the 
sense that they are compatible with the discrete time slicing of our causal geometries 
\29\ I47j . The moves are described in detail in the next Section. 

Topology and causality. Because move 3 is self-dual, the CDT program is using a 
set of seven local moves, being a combination of Pachner moves. By construction, these 
moves do not destroy the topology x of the spacetime, and preserve the global 
proper-time foliation. Moves affecting the spatial slices preserve the topology of the 
space. The two features define the entropy factor. In fact, we know neither what class 
of geometries we should consider in quantum gravity nor the proper measure. The CDT 
model is a trial to define these objects. 

Detailed balance condition. The previous two points are independent of the proba- 
bility distribution of configurations specified by the action P[T] = -^e"^^'^ . Given a move, 
it transforms a configuration A into configuration B. In order to ensure that the Monte 
Carlo algorithm probes the configuration space with probability P[T], the move has to 
be accepted with some probability weight — )• B) and satisfy the detailed balance 
condition: 

P{A)W{A -^B) = P{B)W{B A). 

This condition determines the weights, however not uniquely. The probability P{A) de- 
pends on the considered model. The transition probability — )• is a product of two 
factors: 

W{A ^B) = W^^\A ^ B) ■ W^^\A B). 

The first Vl^^^^(^ — t- B) is determined by the way we choose a position in A where the 
move is performed and by the move type. The second W^'^\A — ?• B) is chosen such 
that the detailed balance condition is satisfied. Once we have decided what Monte Carlo 
moves, and how we implement them, the weight iy^^^(^ — t- B) is fixed and depends on the 
internal probability of choosing a move. A good Monte Carlo algorithm has the weights 
VF^^^(^ —7- B) and W^'^\B — t- A) close to 1. This is equivalent to a high acceptance of 
moves. 

Autocorrelation and thermalization time. An accurate approximation of the av- 
erage ()6.2p requires a suitably large sample of independent configurations {T"^^^ ) • • • , 7"^^^}. 
Configurations are collected in intervals, measured in a number of Monte Carlo steps, 
whose length determine correlations between configurations. Intuitively, the larger is the 
number of Monte Carlo moves which separate triangulation T^^'^ and T^^'^'^'^ the more 

^The numeration of moves corresponds to the convention used in the computer program. 
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Figure 6.1: A two-dimensional visualization of a starting configurations. It presents a 
minimal triangulation of a thin tube with glued ends, thus having a topology of a torus 
X S^. In four dimensions the topology of our interests is x S^. 

uncorrelated they are. The separation should be chosen properly based on the autocorre- 
lation time, which can be estimated using some slowly changing variables, like a width of 
the volume distribution N(i) or its mass center position. We will refrain from discussing 
the details of this numerical procedure. Before we can start measurements, the config- 
urations have to undergo a thermalisation procedure, which consists on sufficiently long 
execution of the algorithm without performing measurements. 

Tuning of K4. The Regge action (|1.5p depends on three bare coupling constants 
Kq, a and K4. For technical reasons of simulations it is preferable to enforce that the total 
number of four-simplices N41 fiuctuates around some prescribed value nvolume during the 
Monte Carlo simulation. The coupling constant K4 is related to the bare cosmological 
coupling constant. If the value of K4 is smaller than its critical value K^, the total volume 
N4 tends to explode. In the opposite case K4 > K^, the configuration collapses. Thus 
K4 should be tuned to its critical value. In practice we use a modified Regge action by 
adding a term +e|A'4i — nvolumej which, for K4 w K^, ensures that A'41 fluctuates around 
nvolume. Alternatively, we can add a quadratic term +e(A'4i — nvolume)^ which is also 
studied. Presented results are consistent in both cases. 

Starting configurations. Before we can apply the moves and generate a Markov chain 
we must create the initial configuration. We construct a minimal configuration with a 
topology X S"^. Fig. 16.11 presents a two-dimensional analogue of the starting configura- 
tion. First, we assemble the simplest three-dimensional simplicial complex with a topology 
S^, corresponding to a spatial slice at a fixed time. It is a boundary of a 4-simplex and 
contains five tetrahedra glued to each other. Let us label five vertices of this complex by 
0, 1, 2, 3, 4. Then the five tetrahedral interfaces are denoted as 

[0, 1, 2, 3], [0, 1, 2, 4], [0, 1,3,4], [0, 2, 3, 4] and [1, 2, 3, 4]. 

Next, we take two such structures, first at time i with vertices 0, 1, 2, 3, 4, second at time 
i + 1 with vertices 5, 6, 7, 8, 9. And join them with a minimal number of 4-simplices. Each 
of five interfaces at time i becomes a base of a simplex with a top vertex at time i + 1. 
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thus creating five {4, 1} four-simplices. Let us denote these by 

{0,1, 2, 3, 5}, {0,1, 2, 4, 6}, {0,1,3, 4, 7}, {0,2, 3, 4, 8}, {1,2, 3, 4, 9}. 

These are supplemented by 10 simplices of the type {3, 2}: 

{0, 1, 2, 5:6}, {0, 1, 3, 57f}, {0, 2, 3, M}, {1, 2, 3, 5;9}, {0, 1, 4, 6;?}, 

{0, 2, 4, M}, {1, 2, 4,6:9}, {0, 3, 4,7:8}, {1, 3, 4, 7:9}, {2, 3, 4, 879}. 
Likewise, mirror reflected simplices are created, 10 simplices of the type {2,3}: 

{0, 1, 5:6:7}, {0, 2, 5:6:8}, {O, 3, 57^:8}, {0,4,67^}, {1,2, 5:679}, 
{1, 3, 57^9}, {1, 4, 6:7:9}, {2, 3, 5X9}, {2, 4, eTM}, {3, 4, 7X9}, 
and finally 5 simplices of the type {1,4}: 

{0, 5X7:8}, {1, 5X7:9}, {2, KXM}, {3, KXM}, {4,677X9}- 

The result of these operations is a four-dimensional segment with a topology [0, 1] x 5^. 
In total, one segment is composed of five simplices of type {4,1}, ten simplices of type 
{3,2}, ten simplices of type {2,3} and five simplices of type {1,4} summing up to 30 
simplices. The boundary of a segment consists of two identical components, built of five 
interfaces. A number of segments can be iteratively joined one after another forming a 
simplicial manifold of a desired time length. In order to obtain a starting configuration 
with a topology x S^, i.e. periodic boundary conditions, the interfaces at the beginning 
and at the end are glued together. 

Algorithm 

The scheme of the algorithm is as follows: 

1. We select the type of a move at random. The probabilities are chosen in such a way 
that the numbers of accepted moves of each type are approximately the same. 

2. Each move type is related to a characteristic structure of simplices, which must be 
met in order to be accepted. We choose at random a point, link, triangle or simplex 
that satisfy these conditions. The design of the program provides a quick access to 
such elements. 

3. We check if the new simplicial complex corresponds to a triangulation of a manifold. 
One must avoid creating simplices or subsimplices already present in the triangula- 
tion. Also, the move must not destroy the foliation. 

4. We compute the Metropolis weights W^"^^ and accept or reject the move with proper 
probabilities. 

5. We update the lattice structures. 

In the next Section we describe all types of moves, in particular we describe the con- 
ditions which must be satisfied in order to perform a move. 
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6.3 Monte Carlo Moves 

We start the discussion of Monte Carlo moves with recalhng the Regge action which 
determines the probabihty distribution (|6.ip on the ensemble of causal triangulations. The 
implementation of moves may be found in files moves/move?? . c (?? denotes a particular 
move name). The action used in the program is parametrized as: 

S[T] = -{Ko + 6A)A^o + K^N^ + AA^4i, 

where Kq, K4 and A are bare coupling constants and Nq is the number of vertices, 
is the number of simplices, A^4i is the number of simplices of the type {4, 1}. In the 
remaining part of this Section we focus on particular types of moves. 



Move 2 

Move 2 replaces the interface (0,1,2,3) between the two four-simplices by a link (4,5) 
dual to the interface. Here the numbers denote labels of vertices involved in the move, 
and the bracket (...) denotes a tetrahedron built on given vertices. The new link has 
a coordination number 4. We can only create a new time (and not space) link with a 
coordination 4. This can easily be checked before the move is attempted. The first step in 
the move is to pick a {3, 2} or {2, 3} four-simplex from the list and to chose the interface. 
The neighboring four-simplex can be either again a {3, 2} or {2, 3} four-simplex or a {4, 1} 
or {1,4} four-simplex. The move is accepted with a different weight in these two cases, 
since in the first case there are two ways it can be realized and in the second only one way, 
this affects VF^^^ The operation of the move 2 can by schematically written as 

{0,1,2,3,4} + {0,1,2,3,5} 

{0,1, 2, 4, 5} +{0,1, 3, 4, 5} + {0, 2, 3, 4, 5} + {1, 2, 3, 4, 5} 
The effect of this move is: 

• The number of points does not change. AA'^o = 0. 

• New time link (4, 5) with coordination number 4 is created. 

• Four new time triangles with a coordination number 3 are created. 

• Two simplices of type {2, 3} are created. AA''i4 = 0, AA'^23 = 2. 
The weight for the move depends on the Regge action and is given by 

' - PiA) 

where A denotes a triangulation before performing the move and B after performing the 
move. This relation determines the value of W^'^\A —>■ B). 
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Move 2 

It is an inverse of move 2. It requires a (time) link with a coordination number 4. Its 
action can be summarized by 

{0,1, 2, 4, 5} +{0,1, 3, 4, 5} + {0, 2, 3, 4, 5} + {1, 2, 3, 4, 5} 
{0,1,2,3,4} + {0,1,2,3,5} 

• The number of points does not change. AA'^o = 0. 

• Time link [4, 5] with a coordination number 4 is removed. 

• Four time triangles with a coordination number 3 are removed. 

• Two simplices of type {2, 3} are removed. = 0, AN23 = —2. 
Analogically as in previous move, the weight for the move is W2 = e^^". 

Move 3 

A self-dual move. It requires a time triangle (0,1,2) with a coordination number 3, and 
replaces it by a dual triangle (3, 4, 5) with the same coordination number. A triangle 
with a coordination 3 can only be a time triangle. The spatial triangles necessarily have 
coordination 4 or more. It can be easily seen by a following argument: a spatial triangle 
is an interface between two spatial tetrahedra. Therefore it will belong to at least 4 four- 
simplices, for which these tetrahedra are bases and the tips point up or down. In other 
words creating spatial triangles with a coordination 3, although superficially possible, 
creates a topological defect. 

{0,1, 2, 3, 4} + {0,1,2, 3, 5} + {0,1, 2, 4, 5} 

{0, 1, 3, 4, 5} + {0, 2,3, 4, 5} + {1, 2, 3, 4, 5} 

• No new points nor links are created. AA'^o = 0. 

• Triangle (0,1,2) is removed. 
Triangle (3, 4, 5) is created. 

• The number of simplices does not change. = AA^4 = 0. 
All changes are zero, the weight is ■ws = 1. 
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Move 4 

This move adds a new vertex at a given time slice i. It acts on a configuration of {1,4} 
and {4,1} simplices. Let the two simplices be {0,1,2,3,5} and {0,1,2,3,6}, where bars 
indicate the time position of points. We produce a new point {4} between the two sim- 
plices. The move adds 6 new simplices of type {1,4} to our system, and the vertex {4} 
has a coordination number equal 8. 

{0, 1,2, 3,5}+{0, 1,2,3,6} 

{0, 1, 2, 4, 5} + {0, 1, 3, 4, 5}+{0, 2, 3, 4, 5} + {1, 2, 3, 4, 5} + 
{0, 1, 2, 4, 6} + {0, 1, 3, 4, 6}+{0, 2, 3, 4, 6} + {1, 2, 3, 4, 6} 

• One new point {4} with a coordination number 8 is created. 
This property is used as a signal for the inverse move. AA'^o = 1. 

• Four new spatial links with a coordination number 6 are created. 
Two new time links with a coordination number 4 are created. 

• Six new spatial triangles with coordination number 4 are created. 
Eight new time triangles with a coordination number 3 are created. 

• Six new simplices of the type {1,4} are created AA''i4 = 6, AA''23 = 0. 
The weight for this move is u)4 = 



Move 4 

This move is an inverse of move 4. We pick point {4} with a coordination 8 from the 
list of vertex labels. Such a point may be shared only by simplices of type {1,4} and 
the neighborhood satisfies the necessary constraints. In this move the spatial interface 
(0, 1, 2, 3) is created and we have to check if such interface does not already exist. If yes, 
we abort the move. 

{0, 1, 2, 4, 5} + {0, 1, 3, 4, 5}+{0, 2, 3, 4, 5} + {1, 2, 3, 4, 5} + 
{0, 1, 2, 4, 6} + {0, 1, 3, 4, 6}+{0, 2, 3, 4, 6} + {1, 2, 3, 4, 6} 

{0,1, 2,3,5}+{0, 1,2,3,6} 

• Point 4 with a coordination 8 is removed. ANq = — 1 

• Four spatial links with a coordination 6 are removed. 
Two time links with a coordination 4 are removed. 

• Six spatial triangles with a coordination 4 are removed. 
Eight time triangles with a coordination 3 are removed . 
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• Six simplices of the type {1,4} are removed. AA^i4 = —6, AA''23 = 0. 
The weight for this move is 

Move 5 

This move consists on a restructuring of the triangle between four simphces. This move 
requires a spatial triangle (0, 1, 2) with coordination 4. By construction it is shared by four 
simplices of the type {1, 4}. Vertices 0, 1, 2, 3, 4 have the same time position. Vertex 5 is in 
the previous layer and vertex 6 is in the next layer (or reversed). A starting configuration 
must have 4 four-simplices in the following setup: 

{0, 1, 2, 3, 5} + {0, 1, 2, 4, 5}+{0, 1,2,3,6} + {0, 1, 2, 4, 6} 

{0, 1, 3, 4, 5} + {0, 2, 3, 4, 5} + {1, 2, 3, 4, 5}+{0, 1, 3, 4, 6} + {0, 2, 3, 4, 6} + {1, 2, 3, 4, 6} 

• The number of points doesn't change. ANq = 0. 

• One new spatial link (3, 4) with coordination 6 is created. 
This property is used as a signal for the inverse move. 

• Two new time triangles with a coordination 3 are created. 
Three new spatial triangles with a coordination 4 are created. 
Spatial triangle (0, 1, 2) with a coordination 4 is removed. 

• Two new simplices of the type {1,4} are created. = 2, AN23 = 0. 
The weight for the move is = e~2-'^4-2A_ 

Move 5 

This move is an inverse of move 5. We pick a spatial link (3, 4) with a coordination 6. 
Such a link is always shared by 6 simplices of the type {1,4}. 

{0, 1, 3, 4, 5} + {0,2, 3, 4,5} + {1,2, 3, 4, 5}+{0, 1,3, 4, 6} + {0,2, 3, 4, 6} + {1,2, 3, 4, 6} 

{0, 1, 2, 3, 5} + {0, 1, 2, 4, 5}+{0, 1, 2, 3, 6} + {0, 1, 2, 4, 6} 

• The number of points doesn't change. ANq = 0. 

• Spatial link (3, 4) with a coordination 6 is removed. 

• Move creates a spatial triangle (0, 1,2) with coordination 4, unless it already exists. 
If this is the case, the move is aborted. This property is used as a signal for the 
move 5. 
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Two time triangles with a coordination 3 are removed. 
Three spatial triangles with a coordination 4 are removed. 

• Two simplices of the type {1,4} are removed. = —2, = 0. 

The weight for the move is 

At this point we finish the description of the implementation of the Monte Carlo 
algorithm. We presented the basics necessary to create a program generating Monte Carlo 
configurations in the framework of four-dimensional Causal Dynamical Triangulations. A 
more detailed and technical description is beyond the scope of this thesis. We shall omit 
e.g. issues of the random number generator, preliminary acceptance, features specific to 
causal triangulations which allow to reject a move at the initial stage, and other solutions 
affecting the efficiency of the implementation. 

Although we have devoted only one Chapter to the problem of the causal triangulations 
generator, it should be noted that the software development was the major part of the 
work described in this dissertation. 

A copy of the source code together with the user guide and the documentation is 
attached to this dissertation. The author distributes this package on request through the 
e-mail: atgOth . if . u j . edu . pi 



Conclusions 



The aim of this thesis is to present a comprehensive review of recent results obtained 
within the framework of the four-dimensional model of Causal Dynamical Triangulations 
(CDT). We provided answers to questions like: how does a background geometry emerge 
dynamically, what does it correspond to and how to describe quantum fluctuations around 
the average geometry. 

The model of Causal Dynamical Triangulations is a non-perturbative and background 
independent approach to quantum gravity. The foundations of this model are very simple. 
It is a mundane lattice field theory with a piecewise linear manifold serving as a regular- 
ization of general relativity. The introduction of Wick rotation allows to use very powerful 
Monte Carlo techniques and calculate quantum expectation values of observables. 

Based on the Monte Carlo measurements the CDT model predicts the existence of 
three phases. Introduction of a qualitative order parameter describing typical geometries 
appearing in a given phase, allows to make a one-to-one correspondence with three phases 
of an anisotropic theory of a Lifshitz scalar, being an inspiration for Hof ava-Lifshitz gravity 
model. This two models, in addition to formal similarities involving the imposition of a 
global time- foliation, exhibit the same scale dependence of the spectral dimension. From 
a physical viewpoint, particularly interesting is the phase in which the scale factor as a 
function of time behaves as a bell-shaped distribution spontaneously breaking the time- 
translational symmetry. This phase is called de Sitter phase. 

In order to calculate a meaningful average of the three-volume we took into account 
the nonuniform volume distribution and introduced the procedure which breaks the time- 
translational freedom. This step was crucial to evaluate the average scale factor distri- 
bution and confirmed that a background geometry emerges dynamically. The average 
geometry behaves like a well-defined four-dimensional manifold, both the Hausdorff di- 
mension and spectral dimension equal 4 at large scales. We showed that the background 
geometry corresponds to a Euclidean de Sitter space, i.e. a four-sphere, an isotropic and 
homogeneous solution of the vacuum Einstein field equations with a positive cosmological 
constant. Therefore CDT produces a picture of the Universe with superimposed finite 
quantum fluctuations around the classical trajectory. In terms of the lattice spacing, 
configurations resemble a four-dimension spheroid elongated in the time direction. 

The elimination of a translational mode enables to observe the emergence of the back- 
ground geometry, but as a consequence also allows to investigate quantum fluctuations 
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around the average geometry. Using Monte Carlo methods we measured the covariance 
matrix of scale factor fluctuations. Applying the semiclassical expansion we reconstructed 
the effective discrete action describing the three-volume. In the CDT model no reduc- 
tion of degrees of freedom is introduced. The three- volume distribution is obtained by 
integrating out all degrees of freedom except the scale factor. The semiclassical expan- 
sion takes into account both the path integral measure and the bare action, making the 
analysis truly non-perturbative. The extracted effective action may be identified with the 
discretization of the minisuperspace action with an opposite sign. The effect of the re- 
versed sign is a consequence of the entropy factor. Due to the identification, the effective 
coupling constant can be related to the physical gravitational constant, giving a recipe of 
how to obtain a meaningful continuum limit and expressing the lattice constant in terms 
of physical units. 

Understanding of the geometric nature of three-dimensional spatial slices, interpreted 
as leaves of the imposed global proper-time foliation, may be of crucial importance also for 
other attempts to quantize gravity. Although the measurements show that the Hausdorff 
dimension is equal 3, the measured spectral dimension is only half of this value, suggesting 
a fractal structure of spatial slices. Indeed, the fractality is confirmed by a direct analysis 
of tree structures defined in terms of so-called minimal necks. 

Outlook. There are many unresolved questions. The author would like to list few 
points he finds most interesting: 

• What is the role of simplices of type {3, 2} and {4, 1} in better understanding of the 
critical behavior when approaching the phase transition? 

• What is the order of B-C phase transition and how to obtain a proper continuum 
limit? 

• Whether the Hof ava-Lifshitz scenario is realized in our model, i.e. if the anisotropic 
space-time scaling in the UV-limit is observed? 

• How to improve the Monte Carlo algorithm to increase its efficiency near the phase 
transition? 

A physically relevant theory must describe matter fields. The work on incorporating mul- 
tiple massive scalar fields and massive point particle is already in progress. We developed 
an efficient method of calculating field propagator on a random lattice based on a diffusion 
process. An important question is 

• How to observe gravitational interaction between matter in our model and how to 
establish the relation between the effective coupling constant and the conventional 
Newton's gravitational constant? 



Appendix A 

Derivation of the Regge action 



In this Appendix we derive the Euchdean version of the Regge action which, together 
with the introduced measure, defines the partition function (jl.9p of the model of Causal 
Dynamical Triangulations. For simplicity, we start directly with the Einstein-Hilbert 
action (jl.ip after performing a Wick rotation of the spacetime to the Euclidean region 
with the metric signature {+,+,+,+}. The simplices building piecewise linear manifolds 
are embedded in the Euclidean space M^. Simplices are illustrated in Fig. 11.21 At the end 
we show that the derived action is indeed obtained by a Wick rotation of the Lorentzian 
Regge action. The Euclidean Einstein-Hilbert action is given by \37\ [38] 

S§WM ^ d^x^iR - 2A) = - ^ d'x^ Qfci? - a) . (A.l) 

where k = A = The Regge action is obtained by a direct evaluation of the 

Einstein-Hilbert action at a piecewise linear triangulation. It is straightforward to calculate 
that part of the action which depends on the cosmological constant, as it is proportional 
to the total four-volume. It is given by a sum of contributions coming from all individual 
simplices. We distinguish two types of simplices, namely {4,1} and {3,2}. Properties of 
all simplices within one type are identical and their geometry is completely determined by 
the lengths and at- The total volume is therefore given by 

/ d^x^X = A (iV4iVol^^'i> + Af32Vol^3,2}\ ^ ^^ 2) 

where A^4i is the number of simplices of the type {4, 1}, which have a volume Vol''-'^'^-'', and 
A''32 is the number of simplices of the type {3, 2}, which have a volume Vol^^'^^'. Integration 
of the scalar curvature requires some attention. The curvature is singular and localized at 
hinges (triangles). It occurs that, in the case of triangulations, it can be simply expressed 
in terms of deficit angles around those hinges. Let us illustrate this on a two-dimensional 
example. In this case, the scalar curvature is a distribution with a support on vertices, 
and the angle deficit is calculated around them. A triangulation consisting of 6 equilateral 
triangles sharing one vertex and periodically connected along sides may be drawn on a 
plane R^. Because the plane is flat, so is the vicinity of the common point. The sum 
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Figure A.l: An example of a positive (top) and negative (bottom) deficit angle. The first 
triangulation corresponds to a cone with a positive singular curvature at the peak, while 
the second to a saddle. 

of angles adjoint to this point equals 9 = 2tt, the angle deficit 5 = 27r — 9 = and 
scalar curvature vanishes. If such triangulation is built of smaller number of triangles, as 
shown an the top of Fig. lA.ll it forms a cone and the angle deficit is positive, as is the 
curvature. Likewise, a larger number of triangles gathered around one vertex produces a 
saddle, illustrated an the bottom of Fig. lA.H with a negative curvature. 

A generalization to higher dimensions was proposed by Regge in |43j . In this case 
besides the deficit angles the integral of the curvature depends on a volume of homogeneous 
hinges. It was shown that the integral over a neighborhood O containing one hinge A is 
given by 

-h [ d^x^ R = -k{27r-9A)yolA, 
^ Jo 

where 9/\ denotes the sum of dihedral angles around the triangle A and VoIa denotes its 
area. The integral over the whole manifold A4 is obtained by summing contributions from 
all hinges. By a (p, 3 — p)-triangle we will denote a triangle with p vertices in one slice 
and 3 — p in next slice. A {4, l}-simplex contains two types of triangles: 4 (3, 0)-triangles 
and 6 (2, l)-triangles. A {3, 2}-simplex contains 1 (3, 0)-triangle, 6 (2, l)-triangles and 3 
(1, 2)-triangles. Each simplex shares 10 triangles and contributes to the total dihedral 
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angle 9/\. Using these properties we can write the global integral as 



]-k [ d^x^i2 = -27rA:y Vo1a+ 



A 

feiV4i (4eg'(;|voi 



(3,o) + 6e|^Woi 



+ kN- 



32 



0(3,0)^01(3,0) 



'(2,1) 

+ eegf Jvoi 



(2,1) 



+ 



(2,1) + 3eg'|>Vol(i,2)^ 



(A.3) 



where denotes the dihedral angle between two faces of {q, 5 — (7}-simplex sharing 

a {p,3 — p)-triangle, Vol(p 3_p) is the triangle volume. As before all of these quantities are 
functions of and oj. The sum over all triangles A equals ^°1a = ^(3,o)Vol(3 q) + 
(-^(2,1) + -^(i,2))Vo1(2,i) and -/V^^ 3_p) denotes the number of triangles of appropriate type. 
In order to evaluate the Euclidean Einstein-Hilbert action for a specific triangulation, one 
must know volumes and angles associated with simplices. Let us explicitly write down 
the positions of vertices of the two types of simplices embedded in the Euclidean space 
R^. We incorporate the length of spatial links into the coupling constant and rescale 
the edge lengths so that = 1, Of = ^/a. The fourth coordinate defines the original time 
direction. The vertex coordinates of a {4, l}-type simplex are: 



V5 



1 1 ]_ 

2' 2^3' 2^6' 

1 1 ]_ 

2' 2^3' 2V6' 

J 1_ 

V3' 2V6' 
3 



0, 

{ 0, 
0, 



0, 



8a-3 



The distances between vertices {1, . . . , 4} are equal = 1, and the lengths of links ending 
in vertex 5 are equal at = y/a. Using these values we can calculate the simplex volume 
and the dihedral angle between two three-dimensional hyperplanes passing through a given 
triangle and respectively the two remaining vertices of the simplex 



Vol^^'i^ 

cosieg'oj] 
cos[ef2^'j^|] 



a 



96 
1 



2a - 1 
6a - 2' 



For a non-degenerate simplex we get the condition a > |. Similarly, we find the vertex 
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coordinates of a {3, 2}-type simplex. Now two vertices lie in the upper plane 



^^3 



0, ^, 0, oj 

{ 0, 0,-1^}, 
{ 0, 0, 



Again the links lying in the same horizontal plane have the length = 1, and the rest 
have at = y/a. The volume of a {3, 2}-simplex and corresponding dihedral angles, are now 
given by 



96 

{3,2}, 6« - 5 VS 

rr>{3.2}i 1 1 \/4a - 1 
cosB^-is" = — , , Voif2 i) = , 

r^{3,2}i 4a - 3 V4a - 1 

For a non-degenerate simplex we get the condition a > j^. The area of triangles is the 
same for both types of simplices and depends only on and Of. 

To simplify equations (|A.2|) and (|A.3|) we use relations, specific for Causal Dynamical 
Triangulations, between numbers of appropriate sub-simplices. Despite the fact that the 
simplices are a subset of Euclidean space, the triangulation is equipped with the causal 
structure which is a consequence of the original Lorentzian spacetime and each slice is 
a triangulation of a three-sphere. Because each spatial tetrahedron is shared by exactly 
two {4, l}-simplices, and those tetrahedra itself build a simplicial manifold one gets the 
equality between the number of spatial triangles and {4, l}-simplices 

iV(3,o)=iV^'''^. 

Using the above equality, the Euler identity and the Dehn-Sommerville conditions for a 
general four-dimensional triangulation \29\ [66] we get 

^(2,1) + ^^(1,2) = 2iVo + iV4 + A^32. 

By inserting this relation into the sum of ()A.2p and ()A.3p . we get the Regge action for 
a triangulation. It is given by a very simple linear function of three independent global 
parameters, the number of vertices Nq, the number of simplices A^4 and the number A'41 
of {4, l}-simplices, 

S = -iKo + 6A)No + KiN^ + ANii, (A.4) 
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where the coefficients, called bare coupling constants, are nontrivial functions of the cou- 
pling constants present in the continuum Einstein-Hilbert action 

+ 6A = V4a - lirk, (A.5) 

Vl2^^^ ^ V3, 6a -5 

Ka = A H — —k arccos h 

96 4 6a -2 

1 3 4a 



+ \/4a — Ik ( - arccos — , + - arccos — vr | , (A.6) 

\2 2^/3^^ V4^^ 4 8a -4 ^ ^ 



^ V8a - 3 , r-r 1 Tr\ 

Ka + a = A + vofc arccos — , 

96 V 2^6^^ 2 J 



-2 + 

+ V4^t(5arccos^-|). {A.7) 

Equivalently, one can start with the Einstein-Hilbert action (jl.ip and evaluate it at a 
Lorentzian piecewise linear manifold made up of simplices being a subset of Minkowski 
spacetime, as described in [29] (here we use opposite sign convention = a ■ a1). In the 
Lorentzian case a is negative and the simplex volumes are given by 

VolLor(a) = \/^Vol £;uc(a) 

The same relation holds for volumes of time-like triangles. The equations for angles are 
identical in both cases, in the Lorentzian case angles may take complex values. Inserting 
the Lorentzian counterparts of volumes and angles into the Einstein-Hilbert action and 
taking advantage of equations ()A.2p and (|A.2p . we get the Lorentzian Regge action. The 
further step is to analytically continue a from negative values to positive values through 
the lower-half complex plane (-y/— 1 = —i). For a > ^ the Lorentzian action becomes 
purely imaginary. The Euclidean action defined by relation —S^^'^ = iS^°^ coincides with 
equation ()A.4p . 



Appendix B 

Constrained propagator 



Let us consider a multivariate normal distribution with a T component random vector 
X and zero mean. We denote the full propagator inverse by P. Further we impose a 
constraint on the random vector, i.e. we measure x only when it is orthogonal to q 



P[x] 



1 



m 



_ 1 

e 2 



x^Px 



(5(q^x), q^q=l. 



For definiteness we normalized the eigenvector q corresponding to the zero mode The 
partition function of this simple statistical model with a source term k'^x is given by 

Z[k] = [cf xe-^^^P^+l^^- 5(q^x). 



The propagator of the constrained model is given by a derivative of the constrained par- 
tition function, 

In Z[k] 



{XiXj)c 



dkidkj 



k=o 



In order to derive a relation between the constrained propagator C and full propagator 
C = P~^, we use the integral representation of the delta function. The function Z[k] is 
given by 

Z[k] = [ d^xdA e-i-^P-+k^-+^^q^- 



detP 



g|(k+iAq)-C(k+aq) 



'{27r)T / 27r 



det P V q^Cq 



exp 



-k^p-^k ^k^C qcfC k 

2 2q-Cq 



and the constrained propagator C as a function of the full propagator C is given by 

1 



C = C 



TV AC 



CAC, A = qq^. 



(B.l) 
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The matrix P is defined as an inverse of C on a subspace ortliogonal to q, 

PC = 1-A. (B.2) 

It is easy to clieck tliat 

P = (1- A)P(1- A), (B.3) 

solves equation ()B.2p . where we used relation (jB.ip . The constraint inverse propagator P 
can be expressed by the constrained propagator in the following way 

P = (C + A)-^-A. (B.4) 

From (jB.ip and ()B.3p it follows that P and C act in a space orthogonal to q 

Pq = Cq = 0. 
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